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CHAPTER  I 


INTRODUCTION 

In  1956  Townes  et  al  /  proposed  and  successfully  produced 
an  oscillator  which  utilized  the  induced  emission  from  an 
excited  quantum  energy  state.  It  was  named  the  maser.  A 
brief  sketch  of  its  principles  is  as  follows. 

Two  quantum  energy  states  exist.  The  difference  in  their 
energies  is  AE.  In  the  presence  of  an  electromagnetic  field,  a 
system  in  the  lower  energy  state  will  absorb  energy  from  the  field 
at  a  rate  A  which  is  sharply  peaked  when  the  incident  field  is  of  the 
frequency  f  =  ,  where  h  is  Planck's  constant.  This  rate  is  propor¬ 

tional  to  the  square  of  either  the  electric  or  magnetic  field  strength. 
Similarly,  a  system  in  the  upper  energy  state  emits  energy  to 
the  field  at  the  same  rate  A  in  the  same  frequency  and  direc¬ 
tion  as  the  incident  field.  If  an  ensemble  is  composed  of  a 
mixture  of  systems  in  the  lower  and  upper  energy  states,  then 
the  ensemble's  absorption  is  A{nj-n2),  where  n^  is  the  number  of 
particles  in  the  lower  state  and  n^  is  the  number  of  particles 
in  the  upper  energy  state. 

If  we  wish  to  make  an  amplifier,  we  must  obtain  from 
such  a  device  a  field  which  is  larger  than  the  field  which  we  feed 
into  the  device.  In  a  maser,  this  means  that  there  should  be 
a  net  emission  of  energy  ^  the  system  into  the  field.  This 
requires  that  A(nj-n2)  must  be  negative  or  that  n2  >  n^.  If  the 
device  were  lossless,  we  would  have  an  amplifier  whenever  this 
requirement  is  fulfilled.  In  practice,  A(nj-n2)  must  be  sufficiently 
negative  to  overcome  inevitable  losses.  For  an  oscillator, 

A(nj-n2)  must  be  negative  enough  to  overcome  not  only  its  own 
losses  but  the  losses  in  the  load  as  well. 

At  thermal  equilibrium,  the  lower  energy  state  is  always 
more  heavily  populated  than  the  upper  state.  The  ratio  is  given 


-T- 


by  Boltzmann's  distribution 


^  _  g-AE/kT 


Therefore,  some  process  is  required  to  upset  this  thermal  equilibrium 
before  a  maser  can  be  made  to  operate.  This  process  is  called 
inversion  (because  it  inverts  the  population  of  the  states  at 
thermal  equilibrium)  or  pumping. 

In  the  first  maser  due  to  Gordon,  Zeiger,  and  Townes,  a 
beam  of  ammonia  molecules  passed  through  a  multipole  electric 
field  with  a  gradient  which  deflects  high  energy  molecules  different¬ 
ly  from  the  low  energy  molecules.  Thus,  a  separation  of  high  and 

low  energy  molecules  is  achieved.  This  creates  a  beam  of  mole- 
n> 

cules  in  which  >  >  1.  This  beam  satisfies  the  condition  for  a 
“1 

maser.  This  beatn  is  passed  on  into  a  cavity  which  is  tuned  to  the 
transition  frequency  f  =  »  where  h  is  Planck's  constant.  If 

the  Q  is  high  enough  to  give  A  (the  transition  probability)  its 
proper  value,  maser  oscillation  or  amplification  takes  place. 

The  maser  aroused  great  interest  because  of  two  advantages. 

(1)  As  an  oscillator  it  has  very  great  stability. 

(2)  As  an  amplifier  it  has  a  much  lower  noise  figure  than  that  of 
vacuum  tubes  or  transistors. 

These  advantages  were  demonstrated  by  both  theory  and  experi¬ 
ments. 

The  success  of  the  first  maser  naturally  stimulated  vigorous 
efforts  to  improve  its  characteristics.  The  ammonia  maser  had  a 
low  power-handling  capacity.  One  way  to  overcome  this  character¬ 
istic  is  to  use  a  solid  as  the  maser  material.  The  higher  density 
of  active  atoms  would  allow  both  higher  power  capacity  and  com¬ 
pactness.  Zeeman  levels  of  paramagnetic  ions  in  solids  provide 

convenient  energy  levels  for  masers  in  the  microwave  frequency 
4,  5 

range. 
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In  such  a  maser,  a  dielectric  solid,  containing  a  small 
amount  (typically  .  05%)  of  paramagnetic  centers,  is  subject  to 
a  strong  d.  c.  magnetic  field.  The  magnetic  moment  due  to 
electron  spin  interacts  with  the  magnetic  field  to  give  an 
energy  gPlT-  "S-  ,  where  g  and  p  are  constants,  H  is  the  d.  c, 
magnetic  field,  and?  is  the  electron  spin.  Discrete  quantum 
energy  levels  result  from  the  discrete  quantum  orientations 
of  S  with  respect  to  H.  Such  levels  can  be  inverted  by  several 
methods  developed  in  the  study  of  nuclear  magnetism.  Examples 
are  adiabatic  fast  passage  (AFP),  180°  pulse  inversion  or  d.  c. 
field  reversal. 

These  masers  can  generate  a  moderate  amount  of  peak 
power  but  have  the  serious  drawback  that  they  can  only  be 
operated  a.s  pulsed  devices.  This  is  because  the  inversion  process 
and  maser  process  can  never  be  simultaneous  .  First,  an  inver¬ 
sion  process  is  applied.  When  this  process  is  over,  the  maser 
process  takes  over.  The  maser  process  can  last  only  as  long  as 
the  state  of  inversion  exists.  As  the  maser  process  proceeds, 
inversion  is  reduced  until  the  maser  process  is  no  longer  possible. 
The  system  must  then  be  allowed  to  relax  to  thermal  equilibrium 
before  the  inversion  process  is  applied  again.  Because  of  this 
limitation,  the  two-level  solid  state  maser  is  considered  now 

only  as  a  measure  to  reach  the  millimeter  wave  region  by  pulsed 

2  3 

d.  c.  magnetic  fields.  ' 

4 

In  1957  Bloembergen  proposed  a  new  type  of  maser  which 
could  be  operated  continuously.  Briefly,  the  principle  was  as 
follows. 

An  atom  has  three  quantum  energy  states  with  energies 
Ej,  E2,  and  E^,  respectively.  The  energy  differences  are 
AE2^>  ^^32  ^^31  thermal  equilibrium, 

the  probability  of  occupation  of  the  three  levels  is  given  by 
Boltzmann's  distribution.  Each  higher  level  has  a  smaller 
population  then  its  lower  neighbor. 
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Now(  if  a  strong  electromagnetic  field  of  frequency  fg'  *  is  applied 
induced  transition  causes  atoms  at  to  be  excited  to  and  atoms  at 
E^  to  be  de-energized  to  E^.  This  is  called  pumping.  The  induced 
transition  rates  per  atom  for  these  two  processes  are  equal  for 
these  equally  degenerate  states.  The  net  rate  is  then  given  by 


A(n^-n^),  where  A  is  induced  transition  rate  and  n^  and  n^  are  the 
number  of  units  at  energy  E^  and  E^,  respectively.  If  this  rate 
A  is  very  much  larger  than  the  mechanisms  which  promote 


thermal  equilibrium,  a  steady  state  is  achieved  when  A(nj-n2)  =  0. 
Therefore,  with  strong  inducing  fields,  n^  g  n^.  This  condition 
is  called  saturation.  If,  furthermore,  A  is  very  much  larger  than 


the  mechanisms  (relaxation)  which  bring  E^  and  E^  into  thermal 
equilibrium  with  E2,  and  these  mechanisms  are  of  equal  strength 


for  and  E^,  the  population  of  E^  will  remain  approximately 

at  its  value  before  pumping.  This  value  is  n^  =  n°e  ^^2l/kT^ 

where  n,  is  the  thermal  equilibrium  value  of  n. . 

^  ^  „  -AE,,/kT 

n°  +  n°e  n°  -AE^^/kT 

After  pumping|n^  - -  g  -^(l  +  e  ). 


If  AE3j/kT<<l,  n2  =  n°(l- AE2j/kT),  n^  =  n^  g  n°/2  (Z-AEj^/kT) 

=  n” (1  ~  AE^j/ZkT).  Therefore,  inversionis  achieved  between  levels 
E3  and  E2  if  <  2AE2^-  If  >  2AE2^>  inversion  is  achieved 
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between  levels  £2  and  £^.  If  the  relaxations  to  £2  from  £^  and  £ 
are  not  equal,  the  relations  are  modified  to  give 


"“23^^32  “21^^21 


"3  ■  "2  =  -wr~ '  ;:i23+‘^2  + 


where  u..  is  the  relaxation  rate  from  level  i  to  level  j.  Bloembergen 
suggested  5%Ni95%ZnSiF^  •  6H2O  and  l%Gd991tLa  {C2H5SO^)3- 9H2O 
as  possible  material  for  hie  maser. 


The  three-level  maser,  as  Bloembergen's  device  is  called, 

was  first  successfully  operated  by  McWhorter  and  Myers^  using 

3+  7 

the  Zeeman  levels  of  Cr  in  KCr(CNz).  Kikuchi  et  al,  constructed 
a  more  successful  model  using  the  Zeeman  levels  of  Cr  in 
AI2O3  (ruby).  Ruby  masers  are  now  in  practical  operation  as  radio- 
astronomical  preamplifiers  at  microwave  frequencies. 

The  success  of  the  ruby  maser  naturally  turned  the  atten¬ 


tion  of  researchers  to  other  materials  which  might  be  suited  for 
three-level  masers.  Ruby  is  an  excellent  material.  However, 
practical  operation  is  limited  to  below  approximately  10  kmc. 
There  exists  a  need  for  oscillators  and  amplifiers  at  higher  fre¬ 
quencies.  It  is  also  desirable  to  operate  masers  without  liquid 

>|c 

helium.  The  search  for  new  materials  is  to  an  extent  guided  by 


Maiman  (70)  in  Quantum  £lectronics  Conference  has  reported 
that  ruby  masers  have  been  successfully  operated  at  liquid  nitro¬ 
gen  and  dry  ice  temperatures.  There  are  no  reports  of  its  prac¬ 
tical  use.  It  should  be  noted  that  if  a  maser  material  is  suited  for 
operation  at  liquid  nitrogen  temperature,  then  it  can  be  made  to 
perform  even  better  at  liquid  helium.  Therefore,  it  is  possible 
that  masers  will  always  be  operated  at  liquid  helium  temperature 
as  long  as  it  is  available.  The  cryogenics  involved  is  expensive 
and  inconvenient  from  the  point  of  view  of  technical  applica-tion. 
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these  two  considerations.  So  far,  the  search  has  turned  up  five 

8  9  10  11 

pronnising  materials.  They  are  beryl,  TiO,  (rutile),^'  ’ 

^  12  13 

MgWO^,  ZnWO^,  and  CaWO^  (calcium  tungstate  or  scheelite  ’  ). 

The  first  two  materials  are  capable  of  operation  well  into  the 
millimeter  wave  region.  The  tungstates  are  suitable  for  use  in  the 
10  -40kmc  range.  Sn02  is  also  reported  to  be  under  investigation. 

CHAPTER  II 

PARAMAGNETIC  RESONANCE  IN  CRYSTALS: 

A  REVIEW 

The  susceptibility  of  solids  differs  from  that  of  free  atoms 
(gas)  in  that  the  atoms  are  ionized  in  the  solid  and  that  these  ions 
have  considerably  modified  orbital  states  as  the  result  of  the  elec¬ 
tric  field  set  up  by  the  other  ions  in  the  crystal.  Most  ions  have 
complete  shells.  Thgis,  they  are  diamagnetic.  The  ions  with 
incomplete  shells  can  be  paramagnetic.  If  they  were  not  subject 
to  the  crystalline  electric  field,  the  paramagnetism  would  have 
contributions  from  both  the  orbital  and  spin  magnetic  moments 
of  the  ion.  However,  in  many  cases  the  crystalline  field  "quenches" 
the  orbital  moment  and  leaves  a  paramagnetism  which  is  due  to 
spin  alone. 

The  energy  levels  used  in  the  three-level  maser  are  the 
Zeeman  levels  of  the  lowest  orbital  state  of  a  paramagnetic  ion 
in  a  solid.  (Higher  energy  levels  cannot  be  used  because  only  a 
very  small  percentage  of  the  ions  are  in  such  states,  particularly 
at  low  temperatures.  )  The  study  of  maser  materials  is,  therefore, 
the  study  of  the  lowest  energy  levels  of  paramagnetic  ions  in 
crystals. 

The  theory  of  paramagnetism  in  solids  was  pioneered 
by  Van  Vleck^^  in  the  1930's.  Considerable  theoretical  develop¬ 
ment  also  took  place  in  Holland  in  the  1930's  and  1940's.  In  fact, 
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a  firm  basis  of  the  present  theory  was  well  established  during  this 
early  period.  It  is  a  very  remarkable  achievement  when  one  con¬ 
siders  the  fact  that  the  only  source  of  experimental  data  were  sus¬ 
ceptibility  measurements.  They  were  not  able  to  observe  the  energy 
levels  directly  by  microwave  resonance  as  we  do  today. 

A  new  approach  to  the  problem  was  opened  by  Zaviosky^^ 
in  194^  when  he  was  able  to  observe  microwave  paramagnetic 
resonance  of  electrons.  Since  then,  the  field  has  progressed  very 
rapidly.  The  major  part  of  the  work  in  the  field  in  the  last  ten 
years  has  been  in  England,  Good,  organized  accounts  of  the 

17 

present  approach  to  the  theory  are  given  by  Bleaney  and  Stevens 
18 

and  by  Low.  Experimental  results  have  been  collected  by  Bowers 

20  21  19 

and  Owens  and  by  Orton.  More  recently,  Koster  has  developed 

a  new  theory  based  on  group  theory.  The  general  theory  of  atomic 

22 

spectra  is  covered  completely  by  Condon  and  Shortley,  The  theory 

of  transitions  metal  spectra,  in  particular,  is  covered  by  J.  S. 

78 

Griffith.  These  theories  are  retraced  in  some  detail  in  this 
chapter. 


II.  1  The  Hamiltonian  for  a  Paramagnetic  Ion  in  a  Crystal 


Any  discussion  of  quantum  energy  levels  must  begin  with  the 
Hamiltonian.  The  Hamiltonian  for  a  paramagnetic  ion  in  a  c  ry stal  is  given 
by  V-  P“  r,.Z  V-  .2 


+  ^cryst.  + 

i  ^i  i  >  j  ij  ij 


Y  \..i.s. 

Li  ij  1  J 


t  Jlf- U.*zs. 


(2.1) 


B 
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The  first  term  is  the  kinetic  energy  of  the  electrons  and  the 
potential  energy  of  the  electrons  in  the  field  of  the  nucleus.  The 
second  term  is  the  potential  energy  of  the  electrons  in  the  fields  of 
each  other.  The  third  term  is  the  potential  energy  of  the  electrons 
in  the  field  of  the  surrounding  crystal.  The  fourth  is  the  spin-orbit 
coupling  energy.  The  last  term  is  the  Zeeman  energy. 

This  Hamiltonian  is  too  complex  to  be  solved  exactly  for  the 
eigenstates  and  eigenenergies.  In  order  to  use  perturbation  methods, 
the  relative  magnitudes  of  the  various  terms  in  the  Hamiltonian 
must  be  established.  The  largest  terms  are  the  first  and  second 
terms.  The  relative  magnitudes  of  the  remaining  terms  depend  on 
the  ion  in  question  and  its  environment. 

For  the  rare  earths,  it  is  agreed  that  the  order  of  the 
interactions  in  descending  magnitude  is  electron-electron 
interaction  (called  Coulomb  interaction),  spin-orbit  coupling, 
followed  by  the  crystalline  field  interaction  and  the  Zeeman 
energy.  The  crystalline  field  has  a  relatively  weak  effect  because 
the  electrons  which  give  rise  to  paramagnetism  in  the  rare 
earths  are  shielded  by  an  outer  complete  shell. 

In  the  iron  group  the  situation  is  not  clear  cut.  The 
electrons  which  contribute  to  paramagnetism  are  the  outer¬ 
most  electrons.  They  feel  the  full  effect  of  the  crystalline  field. 
Consequently,  the  crystalline  field  interaction  is  comparable 
to  the  electron- electron  Coulomb  interaction.  In  fact,  in  many 
cases  it  exceeds  an  important  part  of  the  Coulomb  interaction  so 

that  researchers  have  differed  in  their  approach  to  the  problem. 

17 

Bleaney  and  Stevens  treat  the  crystalline  field  as  a  perturbation 

on  the  Coulomb  interaction.  This  same  approach  was  taken  by  Van 
15 

Vleck£t^.  in  the  early  analysis  of  alums.  This  approach  is 
outlined  in  some  detail  in  this  chapter. 

Other  researchers  have  approached  this  problem  by 
considering  the  crystalline  field  to  be  more  important  than  the 
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non-central  field  part  of  the  Coulomb  interaction.  This  is  some¬ 
times  referred  to  as  the  crystalline  field  representation.  Examples 
are  the  presentation  of  paramagnetic  resonance  by  Bowers  and 
Owens.  In  ruby,  for  example,  the  optical  spectrum  is  discussed 
more  readily  by  the  crystalline  field  representation.  Practically 

all  current  literature  on  the  optical  spectrum  of  ruby  is  given  in 
23 

this  manner.  Low  has  also  discussed  this  problem.  He  finds 
that  the  crystalline  field  is  of  such  magnitude  that  it  causes  a 
considerable  mixing  of  eigenstates  in  the  LS  coupling  scheme  where 
representation  is  diagonal  in  the  Coulomb  interaction  Hamiltonian. 

Since  we  are  interested  mainly  in  the  paramagnetism 
of  the  iron  group,  the  case  of  the  rare  earths  is  omitted  from 
further  discussion.  We  will  now  proceed  with  the  discussion  of 
the  energy  levels  in  the  iron  group. 

Even  with  perturbation  methods,  the  total  Hamiltonian 
(2.  1)  is  still  too  complex  to  solve  exactly.  Therefore,  a 
simplifying  approximation  is  made  in  the  treatment  of  the 
Coulomb  interaction. 

II.  2  The  Configuration 

The  simplification  is  to  replace  the  electron-electron 
electrostatic  energy  as  an  equivalent  central-field  potential 
energy.  This  is  not  an  inverse  square-law  type  potential.  It 
decreases  much  more  rapidly  with  increasing  r  (the  distance 
from  the  nucleus).  When  this  is  done,  the  term  can  be  lumped 
with  the  nucleus-electron  energy  to  give 

V 

<2.2) 

i 

for  the  first  three  terms.  These  two  terms  are  the  largest 
terms  in  the  Hamiltonian.  Therefore,  this  partial  Hamil¬ 
tonian  can  be  solved  for  eigenstates  and  eigenenergies  which, 
in  turn,  can  be  used  as  the  basis  for  perturbation  calculations 
of  the  states  and  energies  for  the  complete  Hamiltonian. 
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This  partial  Hamiltonian  does  not  involve  electron-electron 
interactions.  Thus,  the  equation  is  separable  in  the  coordinates 
of  each  electron.  The  eigenfunction  is  the  product  of  single  elec¬ 
tron  eigenfunctions,  and  the  energy  is  the  sum  of  individual  elec¬ 
tron  energies. 

The  single -electron  eigenfunction  in  the  central -field  poten¬ 
tial  without  an  inverse  square  dependence  on  r  differs  from  the  in¬ 
verse  square  (hydrogen  atom)  case  in  that  the  energy  is  a  function 
of  t  ,  the  orbital  angular  momentum  quantum  number.  The  quan¬ 
tum  numbers  are  still  n  (radial  quantum  number),  i  (orbital 
angular  momentum  quantum  number),  m^  (z-component  of  orbital 

angular  momentum),  and  m  (z-component  of  spin). 

s 

The  group  of  all  states  having  the  same  quantum  numbers  n 
and  I  is  called  a  shell.  A  shell  is  of  particular  significance  be¬ 
cause  of  several  reasons.  When  all  states  in  the  shell  are  occupied, 
the  electrostatic  potential  due  to  the  electrons  in  the  shell  becomes 
isotropic  (spherically  symmetric).  The  total  orbital  and  total 
spin  angular  momenta  of  the  electrons  in  the  shell  become  zero. 
Therefore,  ions  with  full  shells  are  diamagnetic. 

The  eigenfunction  for  atoms  or  ions  with  more  than  one  elec¬ 
tron  can  be  formed  as  a  combination  (product)  of  single  electron 
eigenfunctions.  We  are  interested  in  the  ground  state  (lowest 
energy  eigenfunction).  This  might  be  formed  by  taking  the  lowest 
energy  state  for  each  electron,  were  it  not  for  Pauli's  exclusion 
principle  which  prohibits  such  a  combination. 

The  exclusion  principle  by  Pauli  states  that  no  two 
electrons  can  occupy  identical  eigenstates.  An  alternate  state¬ 
ment  is  that  electron  eigenfunctions  must  be  anti- symmetric  with 
respect  to  permutations  of  the  electrons. 

The  lowest  energy  state  of  a  many-electron  atom  is  then 
constructed  by  filling  the  single -electron  eigenstates  of  the 
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atom  one  by  one,  starting  with  the  lowest  energy 
state  and  moving  upward  in  energy  until  all  the  electrons  are 
assigned  to  a  state.  This  means  that  we  start  with  states  of 
lowest  n  first.  Among  states  of  the  same  n,  the  states  of  lowest 
I  usually  have  the  lowest  energy.  Such  combination  states  are 
called  configurations.  Restated,  a  configuration  describes  the 
state  of  an  atom  by  specifying  the  state  of  each  individual  electron 
without  specifying  any  relation  between  them.  The  energy 
difference  between  configurations  is  usually  very  large.  Therefore, 
all  other  terms  in  the  Hamiltonian  are  considered  to  be  effective 
only  within  a  configuration. 

II.  3  The  Term 

An  approximation  of  the  electron-electron  electrostatic 
interaction  by  a  central-field  potential  led  to  a  spherically 
symmetric  central-field  potential  for  a  full  shell.  We  have 
arrived  at  a  consistent  result  for  this  configuration.  However, 
in  the  case  of  an  incomplete  shell,  corrections  must  be 
made  to  account  for  the  non-central-field  nature  of  the 
electron- electron  electrostatic  interaction.  This  is  done  by 
grouping  together  all  the  electrons  in  complete  shells  and 
the  nucleus  into  a  "pseudo-nucleus"  with  a  central  potential 
The  interaction  of  the  remaining  electrons  (not  in  full  shells) 
with  the  "pseudo-nucleus"  and  with  each  other  is  then  con- 

22 

sidered.  This  is  discussed  in  detail  by  Condon  and  Shortley. 

It  is  too  involved  to  be  repeated  here. 

The  effect  of  the  non-isotropic  electrostatic  field  due 
to  incomplete  shells  is  to  split  the  configuration  into  terms. 

It  is  recalled  that  the  configuration  was  the  classification  of 
an  atomic  state  by  the  description  of  each  electron  in  the  atom. 

This  is  approximately  true.  Deviations  from  this  rule  cause 
the  rare  earths  to  have  incomplete  shells. 


-11  ■ 


It  did  not  specify  the  orientations  of  the  monnents  of  one  electron 

with  respect  to  another.  A  term  is  an  atomic  state  in  which  the 

relative  orientations  of  the  angular  momenta  of  the  electrons  in 

a  configuration  are  specified  to  give  particular  total  moments. 

A  term  is  a  substate  of  a  configuration.  For  example,  the  con- 
2 

figuration  3d  (two  3d  electrons)  can  give  rise  to  a  total  orbital 
moment  which  is  any  integer  between  4  and  0.  The  total  spin 
can  be  either  1  or  0.  Specification  of  a  particular  possible 
value  for  each  of  the  two  momenta  identifies  a  term.  The  lowest 
energy  term  is  the  term  in  which  the  total  spin  is  the  maximum 
allowed  by  Pauli's  principle.  Among  terms  with  this  maximum 
spin,  the  term  with  the  largest  total  orbital  moment  consistent 
with  Pauli's  principle  has  the  lowest  energy.  These  results 
are  also  known  as  Hund's  rules. 

The  terms  are  approximate  eigenstates  of  the  atom  as 
long  as  the  electron-electron  interaction  is  stronger  than  the 
crystalline  field  or  spin-orbit  coupling.  In  the  free  atom, 
the  generation  of  the  terms  is  called  Russel-Saunders 
coupling. 

The  energy  difference  between  terms  is,  again,  usually 

17 

in  the  optical  region.  Bleaney  and  Stevens  adopt  the  approach 
in  which  this  splitting  is  considered  greater  than  that  due  to  the 
crystalline  field  which  follows.  Mixing  of  terms  due  to  the 
crystalline  field  is  not  considered.  As  we  have  mentioned 
before,  this  approximation  is  not  very  good.  The  alternate 
methods  are  discussed  later  in  this  chapter. 

Following  the  Coulomb  interaction,  further  pertur¬ 
bation  and  removal  of  degeneracy  occur  through  the  effects 
of  the  crystalline  field  and  spin-orbit  interaction.  These  two 
effects  are  not  conveniently  treated  simultaneously.  Therefore, 
it  is  necessary  to  determine  which  of  the  two  effects  are 
dominant.  It  has  been  determined  that  in  the  iron  group,  where 
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the  incomplete  shell  is  the  outermost  shell  of  the  ion,  the  electrons 
in  the  shell  are  sensitive  to  the  crystalline  field.  The  ionic  eigen¬ 
states  in  the  crystalline  field  are  considered  as  the  basis  upon  which 
the  smaller  spin  orbit  coupling  are  considered  to  act. 

II.  4  The  Effect  of  the  Crystalline  Field  on  the  Orbital  State 

The  crystalline  field  term  in  the  Hamiltonian  is  due 
to  the  interaction  of  the  electrons  in  the  iron  group  ion  with  the 
electrostatic  field  of  its  charged  neighbors.  Strictly  speaking, 
it  is  the  field  due  to  all  the  other  ions  in  the  crystal.  However, 
it  is  often  estimated  by  taking  the  nearest  neighbors  only. 

Because  the  electrostatic  field  in  question  is  the  field 
due  to  charges  other  than  the  field  of  the  central  ion,  the  field 
at  the  central  ion  can  be  represented  as  a  general  solution  to 
Laplace's  equation  up  to  the  radius  of  the  nearest  neighbor. 

+n 

V  =  I  I  <2- 

n  m=-n 

We  are  interested  only  in  that  part  of  V  which  interacts 
with  the  orbital  states  of  the  central  ion.  Therefore,  we  can 
neglect  any  term  in  V  for  which  there  exist  no  matrix  elements 
among  the  orbital  states.  That  is,  if 


for  all  i  and  j,  the  term  in  the  potential  is  ignored. 

The  atomic  orbital  states  to  which  this  perturbing  poten¬ 
tial  is  applied  are  combinations  of  products  of  single  electron 
functions.  Furthermore,  as  we  have  outlined  in  the  previous 
section,  they  are  characterized  by  a  particular  value  of  total 
orbital  moment  L  and  a  particular  S,  of  the  total  spin. 

A  typical  matrix  element  of  V  is  of  the  form 
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(2.4) 


Each  of  the  single  electron  functions  4^^^  are  spherical 
harmonics  multiplied  by  a  radial  function.  For  the  iron  group, 
all  the  electrons  are  in  the  id  shell.  So  we  have  n=2.  The  matrix 
elements  are  proportional  to 

^  I  f(r)|  Y‘'^'Y5jYj'dT  (2.5) 

From  the  properties  of  spherical  harmonics  we  know 
that  this  integral  is  zero  unless  n  <  4  =  2|  and  m'  =  p+m. 
Furthermore,  if  4<'  ^.nd  2ire  of  the  same  orbital  moment  f , 
the  product  4'''!'  i®  even  function.  Therefore,  for  the  product 
4*'  Y^  4*  to  be  non-zero,  Y^  must  likewise  be  even.  This  limits 
the  interest  in  n  to  just  the  even  n's.  For  the  rare  earths,  where 
I  =  3,  we  must  take  all  the  even  n's  up  to  21  =  6.  This  greatly 
limits  the  terms  in  the  crystalline  field  to  be  considered. 

We  can  obtain  further  information  about  the  general 
form  of  the  crystalline  field  by  consulting  the  crystallo-  - 
graphic  data  on  the  crystal  in  question.  These  sources  provide 
all  the  relevant  distances  to  locate  the  position  of  a  parti¬ 
cular  atom  in  the  lattice  and  the  point  symmetry  about  this 
position. 

The  crystalline  field  about  an  ion  must  have  a  form 
which  is  consistent  with  the  point  symmetry  about  the  ion.  If 

a  p-fold  rotational  symmetry  is  present,  the  crystalline  field 
can  only  contain  terms  y|^,  where  m  is  restricted  to  multi¬ 
ples  of  p.  If  a  horizontal  reflection  plane  exists,  the  values 
of  m  must  be  even  for  the  iron  group;  this  is  necessary  because 
I  is  even.  Many  such  relations  serve  to 
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greatly  simplify  the  crystalline  field.  In  fact,  if  the  rotational 
symmetry  about  an  ion  is  sufficiently  high,  the  only  relevant 
portion  of  thecrystalline  field  is  that  of  axial  (cylindrical) 
symmetry.  An  iron  group  ion  in  a  hexagonal  site  would  be  such 
an  example. 

After  the  simplest  relevant  form  of  the  crystalline  field 
is  found,  the  orbital  energies  and  eigenstates  under  its  influence 
are  calculated.  Given  the  matrix  elements  of  V,  this  can  be 
carried  exactly  for  many  cases.  In  the  cases  of  very  low  symmetry, 
the  use  of  perturbation  methods  may  be  necessary. 

The  explicit  calculation  of  the  matrix  elements  can  be  made 
from  (2.  4).  This  is  a  fairly  tedious  process  because  Ip  are  anti¬ 
symmetric  functions  of  many  electrons.  A  method  which  greatly 

simplifies  the  calculation  of  the  non-radial  part  has  been  pro- 
24 

vided  by  Stevens.  In  this  method,  the  crystalline  field  is 

replaced  by  an  equivalent  operator,  which  operates  on  the  atomic 

orbital  states  |  L,  M  >  instead  of  the  individual  electron  orbital 

states.  It  is  based  on  the  proportionality  which  exists  between 

the  matrix  elements  of  variables  x,  y  and  z  as  vector  operators 

and  the  matrix  elements  of  L^,  L^,  and  in  a  manifold  of 

constant  L.  Tables  giving  the  results  of  explicit  calculations 

25 

using  this  method  are  provided  by  Elliot  and  Stevens  for  all 
the  ions  in  the  iron  and  rare-earth  groups.  The  results  are  for 
spherical  harmonics  of  unit  amplitude.  The  actual  field  ampli¬ 
tude  is  required  for  the  final  evaluation  of  the  matrix  elements. 

In  addition,  it  must  be  multiplied  by  the  radial  part  of  the  matrix 
element. 


where  f(r)  is  the  radial  part  of  the  electronic  orbital  function. 
This  is  a  factor  which  is  not  known  to  any  accuracy.  It  is 
because  of  the  unreliable  nature  of  the  values  for  both  the 
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crystalline  field  strength  and  the  function  f(r),  that  the  energy 
levels  in  the  crystalline  field  are  measured  by  optical  absorption 
instead  of  being  calculated.  However,  the  calculations  do  provide 
values  of  relative  energies  which  are  sufficiently  reliable  to  deter¬ 
mine  which  state  is  the  ground  state  and  the  degeneracies  of  each 
level. 

II.  5  The  Effect  of  Spin-Orbit  and  Zeeman  Energies,  The  Case 
of  the  Non-Degenerate  Ground  Orbital  State 

The  final  step  in  the  determination  of  the  lowest  energy 

levels  with  the  total  Hamiltonian  is  the  application  of  the  spin- 

orbit  coupling  term  and  the  Zeeman  term  in  the  Hamiltonian.  The 

26 

current  approach  to  this  problem  is  due  to  Abragam  and  Pryce, 
who  separated  the£:  piroblem  into  the  case  of  the  non-degenerate 
orbital  ground  state  and  the  case  where  the  orbital  ground  state  is 
degenerate.  We  arrive  at  such  ground  states  from  the  splitting  of 
the  term  by  the  crystalline  field  as  described  in  the  preceding 
section. 

In  this  case,  the  degeneracy  which  exists  before  the  appli¬ 
cation  of  spin-orbit  coupling  and  the  Zeeman  energy  is  just  the 
degeneracy  due  to  spin.  This  degeneracy  is  partially  or  com¬ 
pletely  removed  by  spin-orbit  coupling  and  completely  removed 
by  the  Zeeman  energy,  in  any  case.  Abragam  and  Pryce  have 
developed  a  spin-Hamiltonian  which  is  a  pure  spin  operator 
which  is  applied  to  the  spin  states.  In  such  an  operator,  the 
orbital  contribution  through  spin-orbit  coupling  and  orbital 
Zeeman  energy  are  contained  in  constant  coefficients  which 
can  be  evaluated  by  second-order  perturbation  because  the 
orbital  ground  state  is  non-degenerate. 

The  spin-orbit  and  Zeeman  energy  are  given  by  the 
following  Hamiltonian: 
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K=  M-*  S  +  H-  L  +  gpH*  S  =  L  (X.S  +  H)  +  gpH'  S  (2.6) 

If  we  consider  just  the  spin  states  associated  with  the 
ground  orbital  state,  the  matrix  elements  are  of  the  form 

<  0.  m'  ,1^  I  0,  m^  >  ,  where  |0> 
is  the  orbital  ground  state.  For  terms  in  the  Hamiltonian  invol¬ 
ving  Li,  this  means  that  the  matrix  elements  will  contain  a  fac¬ 
tor  <0|l|0>.  For  non-degenerate  orbital  states,  this  is  always 
zero.  Therefore,  the  lowest  spin  states  are  independent  of  the 
effects  of  orbital  moments.  This  is  approximately  true  for  large 
magnetic  fields. 

In  the  range  of  magnetic  fields  used  for  microwave  para¬ 
magnetic  resonance,  this  is  usually  a  poor  approximation. 

The  reason  for  this  is  that  the  terms  in  the  Hamiltonian 
involving  L  do  have  matrix  elements  connecting  the  ground 
orbital  state  to  excited  orbital  states.  This  mixes  a  small 
amount  of  the  excited  state  into  the  ground, state,  so  that  in  the 
second-order  approximation,  the  energy  of  the  spin  states 
becomes  dependent  on  L.  This  dependence  can  be  quite  large 
compared  to  the  Zeemaii  energy  because  the  spin-orbit 
coupling  constant  is  a  fairly  large  constant- 

The  constant  coefficients  in  the  spin -Hamiltonian  are 
related  to  the  spin-orbit  and  orbital  Zeeman  energy  by  the 
second-order  energy  expression 
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The  total  perturbed  energy  is  then 

(2.  7)  +  g^pTI- ■g  (2.8) 

which  is  usually  rewritten 

where  g  and  D  are  tensors.  The  quadratic  term  in  H  is  disregarded 
because  it  is  not  a  function  of  S.  It  cainnot  contribute  to  the  splitting 
of  spin  states.  By  proper  choice  of  axes,  D  and  g  can  usually  be 
diagonalized. 


H  .  =P(g  H  S  +g  H  S  +  g, H  S  ) 
spin  ^®xxx®yyy  ®zzz' 

+  D  +  D  +  D,sf 

•X  -ic  V  V  z  z 


(2.9) 


Terms  of  higher  power  in  S  can  be  obtained  by  going  to- 
expression  for  higher  perturbed  energy.  However,  we  know 
in  advance  that  the  powers  of  H  and  S  must  always  add  up  to 
an  even  number.  This  is  necessary  if  the  spin- Hamiltonian 
is  to  be  invariant  under  time  reversal.^ 

In  principle,  the  spin- Hamiltonian  can  be  derived  and  the 
coefficients  computed  explicitly  if  the  ground  and  excited 
orbital  functions  and  their  energies  were  known.  In  practice, 
we  have  only  qualitative  information  about  the  orbital  levels 
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and  functions.  Therefore,  only  the  form  of  the  Hamiltonian  is 
chosen  to  fit  the  point- symmetry  of  the  ion  site  and  the  coeffi¬ 
cients  determined  by  measurement. 

In  this  spin -Hamiltonian,  the  dipole-dipole  interaction 
between  spins  and  the  interaction  of  any  nuclear  magnetic  mo¬ 
ment  with  the  spin  and  the  orbit  have  been  ignored.  These  fac¬ 
tors  were  taken  into  account  by  Abragam  and  Pryce  to  give  a 
complete  spin -Hamiltonian  of  the  following  form: 


=■5-  D--S  +  pTi-  g-  S  +  ^-  A-  I  +  I  •  P-  I  I 

(2.10) 

where  g,  D  and  A  are  tensors  of  rank  2,  T  is  the  nuclear  spin. 

The  spin- Hamiltonian  appropriate  to  various  common 
symmetry  are  given  below. 


Trijjjonal 

X.  =  4s(S+l) 

=  P  (g,  t  gx(H,S,+  HySyDf  DS^  ti55(F-a) 

•  [sJ-4s(Stl)s2  +4sf)-'^a(S^(S,\s^+(S^  +  S^)Sj 

Orthorhombic 


M  =P(gHS  +  gHS+gHS)  +  DS^+  D,  sf  + 
s  XX  ®yyy®zzz'  xx  yy  zz 


+  D'S^  +  E(S^-S^) 
z  z  X  y 
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Tetragonal 


H  =pfg„H  S  +gX(H  S  +H  S  )lfDS^+^  rs^-|s(S+l)S^+|s^l 

s  Lllzz^  sx  yyj  z  36l.z7  z7zJ 

a,„4  „4  4, 

+  T^S  +S  +S  ) 
o  X  y  z 

Cubic 

'K^=gPff  ■  St  |(s%s^sj) 

One  additional  simplification  can  be  made  in  many  cases 
of  high  symmetry.  This  is  due  to  the  limited  freedom  of  response 
of  small  spin  manifolds  to  their  environment.  A  spin  of  ^  can 
have  only  two  orientations:  up  or  down.  But  the  even  parity  of  all 
spin -Hamiltonian  always  makes  these  two  directions  equivalent. 
The  only  distinction  which  such  a  limited  spin  manifold  can  make 
is  a  difference  in  the  elements  of  the  g -tensor. 

A  spin  of  1  can  distinguish  a  monoclinic  symmetry  but 
no  more.  A  symmetry  which  is  higher  can  be  represented  only 
by  an  axially  symmetric  or  spherically  symmetric  spin -Hamil¬ 
tonian  with  the  exception  of  the  g-tensor.  A  spin  of  3/2  can  dis¬ 
tinguish  a  symmetry  as  high  as  orthorhombic;  higher  symmetries 
can  be  represented  by  an  axial  or  spherical  Hamiltonian  except 
for  the  g-tensor.  A  spin  of  2  can  distinguish  symmetries  up  to 
and  including  trigonal  symmetry.  Finally,  a  spin  of  5/2  is  ca¬ 
pable  of  distinguishing  any  symmetry  representable  by  spheri¬ 
cal  harmonics  of  order  4. 

For  ions  of  iron  group,  it  is  found  that  the  crystalline 
field  about  the  ion  is  almost  always  approximately  cubic.  The 
reason  for  this  lies  in  the  ratio  of  ionic  radii  of  iron  group 
ions  to  common  anion  neighbors,  particularly  oxygen. 

Relatively  small  distortions  from  cubic  symmetry  are  in  the 
form  of  distortions  along  one  or  more  axes  of  the  cube  pro¬ 
ducing  tetragonal  or  orthorhombic  symmetry  or  a  distortion 
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along  a  body  diagonal  to  produce  trigonal  symmetry.  The  spin- 
Hamiltonian  then  takes  a  cubic  form  with  terms  of  lower  symmetry 
added.  When  these  terms  are  of  sufficiently  low  symmetry,  such 
that  they  can  be  expressed  in  spherical  harmonics  of  the  order  2, 
they  take  on  a  relative  magnitude  disproportionately  large  com¬ 
pared  to  the  cubic  part  of  the  spin-Hamiltonian.  This  is  because 
terms  in  the  spin-Hamiltonian  which  are  equivalent  to  spherical 
harmonics  of  order  2  are  the  result  of  second-order  pertur¬ 
bation;  whereas,  the  cubic  part  must  be  the  result  of  fourth 
order  perturbation.  Thus,  a  small  quadratic  term  in  the  crys¬ 
talline  field  can  produce  a  coefficient  in  the  spin-Hamiltonian 
which  is  much  larger  than  that  produced  by  a  large  cubic 
potential.  Thus,  except  for  very  exceptional  cases,  the  quad¬ 
ratic  terms  in  the  spin-Hamiltonian  are,  by  far,  the  most  im¬ 
portant  ones  when  they  exist.  We  can  discuss  the  qualitative 
behavior  of  the  spin  under  the  spin-Hamiltonian  by  just  taking 
the  quadratic  terms  whenever  the  symmetry  is  non-cubic. 

II.  6  Zero-Field  Splitting 

One  of  the  very  important  features  of  the  spin-Hamil¬ 
tonian  is  that  the  degeneracy  of  the  spin  states  is  partially 
removed  even  in  the  absence  of  a  magnetic  field.  Taking  the 
simple,  axially  symmetric  Hamiltonian  as  an  example,  we 
have 


=  "S-  D-  "5  =  DS^ 
s  z 

This  Hamiltonian  splits  the  spin  states  into  levels  given  by 
2 

e  =  Dm  ,  where  m  is  the  z-component  of  the  spin  moment, 
z  z 

This  gives  levels  determined  by  |  m^  |  and  results  in  doubly 

degenerate  levels  (+ni^,  except  for  the  state  m^  =  0.  If 

the  number  of  electrons  is  odd  so  that  the  spin  is  half-integral, 
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all  the  levels  become  doubly  degenerate.  Kramer  has  proved 
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that  this  must  always  occur  for  an  odd  number  of  electrons 
regardless  of  the  symmetry. 

In  an  orthorhombic  symmetry,  for  example,  the  spin- 
Hamiltonian  is 

DS^  +  E(S^-Sy)  =DS^  +  E/2(sJ  +  S^), 

where, 

S ,  =  S  +  iS  and  S  =  S  -  iS  . 

+  X  y  -  X  y 

This  Hamiltonian  mixes  states  m  +  2  with  the  state  m  .  For 

8—  s 

S  =  i/2,  we  obtain  two  levels  both  doubly  degenerate.  The  ma¬ 
trix  of  is 
6 


+3/2 

-1/2 

-3/2 

+1/2 

3/2 

9D/4 

•/TE 

0 

0 

-1/2 

•/IE 

D/4 

0 

0 

-3/2 

0 

- 

9D/4 

^JiE 

+1/2 

0 

0 

-nHE" 

D/4 

The  energies  are  solutions  to  the  secular  equation 
(|D-t  )(^-€)-  3E^  =  0 


(2.  11) 


They  are 


c  = 


- 2 - 


3E 
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The  spin  states  are  obtained  by  substitution  of  back  into  the 
secular  equation.  They  are 


for  the  upper  and  lower  energies,  respectively. 

The  physical  reason  for  this  energy  splitting  in  zero 
field  is  as  follows.  A  stationary  electrostatic  field  locks  the 
electrons  inthe  lowest  energy  position.  If  there  is  only  one  such 
position,  the  electrons  are  stationary  and  lose  their  orbital 
moment.  Spin-orbit  coupling  re-induces  a.  small  amount  of 
orbital  moment  which,  in  the  second- order  perturbation, 
reduces  the  energy  by  the  spin-orbit  energy.  If  the  ion  were 
situated  in  an  isotropic  environment,  the  amount  of  induced 
orbital  moment  and  the  resulting  spin- orbit  energy  would  be 
equal  for  all  orientations  of  spin.  Thus,  we  would  have  no  zero- 
field  splitting.  If  the  environment  is  anisotropic,  the  amount  of 
induced  orbital  moment  can  differ  for  different  orientations 
of  the  spin.  This  results  in  zero-field  splitting.  Because  of  the 
magnitude  of  the  spin-orbit  coupling  constant  (cJlOO  cm.  ^), 
only  a  smal)  difference  in  orbital  moment  can  result  in 
zero-field  splitting  inthe  range  of  a  microwave  quantum. 

II.  7  The  Zeeman  Energy 

The  energy  which  a  magnetic  moment  acquires  in  a 
magnetic  field  is  called  the  Zeeman  energy,  because  of  the 
quantization  of  the  magnetic  moment  in  the  direction  of  the 
applied  field,  2J+1  discrete  levels  result  from  the  applica¬ 
tion  of  the  magnetic  field.  It  is  represented  in  the  spin- 
Hamiltonian  by  the  term. 
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pTT  •  g  -  “5 


If  the  coordinates  are  chosen  properly  with  respect  to 
the  crystalline  field,  the  g-tensor  is  usually  diagonal.  In  this 
case 


p'H-g-'5=8(gHS  +  gHS  +  gHS) 
e  '»x  XX  y  y  y  z  z  z 


For  a  pure  spin  state  g  =  g  =  g  =  2.  003.  Thus,  we  have 

X  y  z 


pH  •  g-  "S  =  gpH-  s- 


In  the  spin-Hamiltonian,  the  small  orbital  component  in 
the  magnetic  moment  causes  g  to  deviate  from  the  pure  spin 
value.  Mathematically,  this  is  due  to  the  cross-terms  between 
H  and  S  in  the  second-order  perturbed  energy  expression  used 
in  the  derivation  of  the  spin-Hamiltonian. 


L.|k>(\S.  +  pH.)| 


2 


which  gives  for  the  coefficient  of  the  terms  pH^S., 


I  <0|Ljk>|2 


2\ 


This  quantity  is  added  to  2.  OOi  to  form  g. 

Physically,  this  quantity  can  be  explained  in  the  fol¬ 
lowing  way.  To  first  approximation,  the  orbital  moment  is 
zero.  The  presence  of  spin  and  the  application  of  a  magnetic 
field  both  serve  to  induce  a  small  amount  of  orbital  moment 
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by  mixing  excited  orbital  states.  The  spin  interacting  with  orbital 
moment  induced  by  itself  gives  zero-field  splitting.  The  spin 
interacting  with  the  orbital  moment  induced  by  the  applied  magne¬ 
tic  field  contributes  to  g.  Similarly,  the  magnetic  field  inter¬ 
acting  with  the  orbital  moment  induced  by  the  spin  contributes 
an  equal  amount  to  g.  In  fact,  the  deviation  of  g  from  2.  OOi  is 
a  direct  measure  of  the  induced  orbital  moment. 

L.  =  .  {X.S+PH) 


The  close  relation  between  zero-field  splitting  and  the 
deviation  of  g  from  2.  003  has  been  pointed  out  by  Bleaney  and 
Stevens.  Namely, 


(g.  -  2.  003) 

- 2 - 


X. 


X  is  of  the  order  of  100  cm  This  shows  that  a  deviation  of  g 
from  2.  003  of  .  01  can  cause  a  zero-field  splitting  up  to  .  5  cm"^. 
The  zero-field  splitting  is  due  to  the  difference  in  D  in  different 
directions.  Therefore,  it  can  be  expected  to  be  somewhat  smaller 
unless  the  are  of  different  signs. 

In  ions  where  the  ground  term  in  the  configuration  has 
no  angular  momentum  the  physical  basis  for  the  spin-Hamil- 
tonian  is  quite  different,  and  these  relations  may  not  hold. 


II.  8  Hyper  fine  Structure 

The  term's-  A-  T  in  the  spin- Hamiltonian  represents 
the  interaction  of  the  nuclear  magnetic  moment  (spin)  of  the 
paramagnetic  ion  with  its  electrons.  This  is  a  static  dipole - 
dipole  interaction  which  is  governed  by  the  relative  orienta¬ 
tions  and  positions  of  the  two  moments. 
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The  effect  is  to  split  each  electronic  spin  level  into  21+1  levels, 
the  spacing  between  the  new  levels  being  proportional  to  the 
value  M  of  the  spin  state.  As  the  result,  each  allowed  spin 
transition,  Am  =+  1>  is  now  split  into  21+1  separate  transitions. 
The  transitions  are  equally  spaced  in  frequency  with 


Af  =  A 


when  A  is  expressed  in  frequency.  This  group  of  21  +  1  lines 
is  called  the  hyperfine  structure  of  a  line.  This  is  illus¬ 
trated  for  the  simple  case  of  "S-  A-l  =  A  -  "5-  T,  S=2  and 
1=1  in  Fig.  2. 

When  A  is  large  compared  to  the  width  of  the  line,  the 
21  +  1  lines  are  distinct.  They  are  easily  indentified  by  their 
equal  intensities  and  equal  spacings.  This  provides  one  way  for 
the  identification  of  paramagnetic  ions  which  possess  known 
nuclear  spin.  When  A  is  small  compared  to  the  line  width, 
it  simply  contributes  to  the  line  width  as  inhomogenous 
broadening. 

Interaction  with  the  nuclear  moments  of  neighboring  nuclei 
is  also  present  when  such  moments  exist.  This  causes  addi¬ 
tional  splitting  of  the  resonance  line.  The  splitting  is  naturally 
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small  and  consists  of  a  great  many  lines.  Tinkham  was 
able  to  resolve  these  lines  in  the  resonance  of  iron  group 
impurities  in  ZnF^-  More  often,  these  lines  are  not  resolved 
and  contribute  to  the  line  width  as  inhomogeneous  broadening. 

II.  9  The  Case  of  the  Degenerate  Orbital  Ground  State 

When  the  ground  orbital  state  is  degenerate  before 
the  consideration  of  spin-orbit  coupling,  the  spin- Hamiltonian, 
which  was  based  on  a  second- order  perturbation  of  the 
orbital  state,  is  no  longer  valid.  Spin- orbit  coupling  will 
completely  mix  the  degerate  orbital  states.  The  qigenstates 
must  be  obtained  by  exact  diagonalization  of  the  spin- orbit 
matrix. 
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Hyperfine  Spectrum, 
Figure  2 
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The  liberal  mixing  of  the  orbital  states  generally 
results  in  a  large  orbital  moment.  This  has  the  following 
effect  of  the  paramagnetic  behavior  of  the  ion; 

(a)  The  value  of  g  deviates  considerably  from  2. 

(b)  The  zero-field  splitting  becomes  very  large, 
usually  beyond  the  microwave  range. 

(c)  The  relaxation  time  T^  becomes  very  short. 

The  last  two  effects  make  this  case  unfit  for  maser  applications. 

There  is>  however,  one  special  case  where,  in  spite 

of  an  orbital  degeneracy,  the  orbital  moment  is  small.  This  is 

the  case  of  Cr^^  in  an  octahedral  (6-coordination  cubic)  field, 

Fe^"*"  in  a  cubic  (8-coordination  cubic)  field,  and  Cu^^  in  an 

octahedral  field.  These  orbital  ground  states  are  doubly 

degenerate  but  have  the  special  property  that  they  are  non- 
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magnetic.  That  is,  the  matrix  of  L  is  identically  zero 
within  the  doublet.  When  the  field  is  not  strictly  cubic,  the 
states  become  weakly  magnetic,  and  the  degeneracy  is 
removed  except  in  the  case  of  trigonal  symmetry.  Cr'*’^ 
in  the  rhombic  field  of  CrSO^  *  was  found  to  have  a  g  of 

1.  95  for  g^  and  1.  99  for  gx.  ® 

In  a  trigonal  symmetry,  the  states  remain  degeneracy 
but  acquire  a  magnetic  moment  in  proportion  to  the  amount 
of  the  trigonal  distortion.  In  this  case,  it  is  speculated  that  a 
distortion  of  the  symmetry  takes  place,  which  removes  the 
degeneracy  and  lowers  the  energy  of  one  of  the  states.  This 
is  called  the  Jahn-Teller  effect.  ^^The  effect  was  originally 
proposed  for  molecular  complexes.  However,  it  has  been 
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applied  to  iron  group  ions  in  hydrated  salts  by  Van  Vleck. 
Bleaney  et  al.  have  observed  a  spectrum  of  Cu^^  at  low 

3|e 

The  6 -coordination  cubic  and  8- coordination  cubic  have  the 
same  form  of  field  when  expressed  in  spherical  harmonics, 
but  the  polarity  is  reversed.  Thus,  the  lowest  levels  in  one 
become  the  highest  levels  in  the  other. 
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temperature  which  could  be  explained  by  the  Jahn-Teller  effect. 

2+ 

The  case  of  Cu  in  the  hydrated  salts  of  trigonal  symmetry 
(fluosilicates,  etc.  )  is  explained  in  the  following  manner.  The  de¬ 
generacy  in  the  trigonal  field  is  removed  by  a  Jahn-Teller  distortion. 
But  the  lowest  orbital  state  has  three  equivalent  orientations  so 
that  three  separate  spectra  are  observed.  At  high  temperature, 

these  orientations  become  mixed,  and  the  spectrum  becomes 

2+  . 

isotropic.  (The  spectrum  of  Cu  in  a  trigonal  field  would  be 
almost  isotropic  if  the  orbital  ground  state  were  non-degenerate.  ) 

II.  10  Other  Methods  for  the  Analysis  of  Paramagnetic  Ions 
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The  spin -Hamiltonian  of  Abragam  and  Pryce  was  based 

on  a  perturbation  calculation  in  which  the  crystalline  field  is 

weaker  than  the  electrostatic  force  between  the  electrons  in  the 

ion  but  stronger  than  the  spin-orbit  coupling.  This  made  the  ground 

term  in  the  configuration  the  starting  point  in  the  calculation. 

The  spin-Hamiltonian  is  explained  in  another  way  by  Bowers 
20 

and  Owens  in  their  review  of  paramagnetic  resonance.  Their 
starting  point  is  the  configuration.  The  crystal  field  acts  on 
each  electron  in  the  configuration,  splitting  the  orbital  degeneracy 
of  the  single  electrons.  If  n  electrons  of  the  same  shell  form  the 
configuration,  the  electrons  are  assigned  to  occupy  the  n 
lowest  energy  orbitals  of  an  isolated  electron  in  the 
crystal  field.  This  forms  the  ground  state.  The  proper 
degeneracy  in  the  orbital  ground  state  is  obtained  this  way. 

The  form  of  the  Hamiltonian  was  obtained  only  from  the  know¬ 
ledge  of  the  ground  orbital  state  degeneracy  and  the  crystal¬ 
line  field  symmetry.  Therefore,  either  method  yields  the  same 
result.  To  the  inexperienced  researcher,  it  is  somewhat 
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surprising  that  the  second  approach,  which  ignores  the  electro¬ 
static  repulsion  between  the  magnetic  electrons,  gives  the 
proper  ground  state. 

In  optical  work,  however,  the  approach  of  Bowers  and 

Owens  is  the  accepted  one.  This  is  known  as  the  crystalline 

field  approximation  because  it  assumes  that  the  crystalline 

field  is  dominant  over  the  electron-electron  electrostatic 

field.  The  reason  for  this  difference  in  approach  is  that,  as 

far  as  the  ground  orbital  state  is  concerned,  the  Russel- 
22 

Saunders  approach  is  good.  However,  the  higher  energy 
levels  of  the  ground  term  in  the  crystalline  field  approach 
the  lowest  levels  of  the  next  highest  term,  and  considerable 
mixing  occurs.  The  Russel-Saunders  coupling  is  no  longer 
accurate.  Thus,  for  optical  work  it  is  a  poor  approximation. 

Some  light  is  shed  on  this  problem  by  the  recent  work 
23 

of  Racah,  Low^t^.  They  have  undertaken  a  project  of 
.finding  the  energy  levels  and  states  of  iron  group  and  rare- 
earth  ions  in  a  cubic  field.  It  is  to  be  extended  to  lower 
symmetries.  They  have  found,  by  use  of  a  computer,  the 
exact  levels  and  states  starting  from  the  LSJM  representation. 
This  is  a  representation  in  which  the  spin-orbit  coupling  is 
diagonal.  Although  the  crystalline  field  is  stronger  than 
the  spin-orbit  coupling,  this  introduces  no  inaccuracy  because 
the  computer  solves  the  problem  exactly  without  resorting 
to  perturbation  methods.  All  the  terms  in  the  configuration 
with  the  ground  S- state  are  used. 

The  result  shows  that  the  ground  orbital  state  is  very 
nearly  the  state  given  by  the  Russel-Saunders  coupling.  Most 
of  the  excited  states  caxmot  be  represented  even  approximately 
by  the  Russel-Saunders  coupling. 
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11.  11  Koster's  Method 


19  34 

Koster  ’  has  developed  a  group-theoretical  method 
by  which  the  number  of  independent  terms  required  in  a  com¬ 
plete  spin- Hamiltonian  is  given. 

In  this  approach,  it  is  assumed  that  the  eigenstates 
and  energies  for  the  complete  Hamiltonian,  in  the  absence  of 
a  magnetic  field,  is  solved.  This  includes  all  the  terms  in  the 
complete  Hamiltonian  including  the  crystalline  field  and 
spin-orbit  coupling.  The  Zeeman  energy  term  is  then  applied 
as  a  perturbation.  This  method  predicted  many  more  inde¬ 
pendent  terms  than  those  included  in  the  spin- Hamiltonian 
by  Abragam  and  Pryce.  For  example,  for  S  =  3/2,  the  spin- 

Hamiltonian  of  Abragam  and  Pryce  terminates  with  terms  in 

2  3 

S.  .  Koster's  method  includes  the  terms  H.S.  . 

1  11 

We  note  that  such  terms  are  allowed  by  the  time- 
reversal  symmetry  of  the  spin- Hamiltonian.  In  fact,  they 
would  have  appeared  in  the  original  Hamiltonian  if  the  per¬ 
turbation  were  carried  out  to  the  fourth  order.  But  so  would 

2  2  3 

a  number  of  other  terms  of  the  form  H  S  and  H  S.  Why 
are  these  terms  not  predicted  in  Koster's  complete  spin- 
Hamiltonian? 

The  reason  lies  in  the  fact  that  Koster's  method  still 
had  to  resort  to  the  perturbation  method  for  the  final  step. 

It  is  a  first-order  perturbation  expression.  Higher  order 
perturbation  is  necessary  to  produce  terms  of  higher  power 
in  H.  Physically,  omission  of  higher  order  perturbation  in 
the  orbital  Zeeman  energy  means  that  the  induction  of 
orbital  moment  by  the  magnetic  field  is  ignored.  This  is 
actually  a  very  good  approximation.  The  spin- orbit  energy, 
which  is  the  next  least  powerful  term  in  the  total  Hamiltonian, 
is  at  least  100  times  stronger  than  the  orbital  Zeeman  term 
at  practical  magnetic  fields. 
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Experimental  evidence  of  the  additional  energy  terms 
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predicted  by  Koster  has  been  reported.  The  effect  is 
small.  It  appears  that  in  most  cases  they  can  be  neglected. 

II,  12  The  Theory  of  T^ 

When  the  spin  system  is  not  in  thermal  equilibrium 

with  a  thermal  bath  with  which  it  is  in  contact,  it  will  relax 

into  equilibrium  in  a  characteristic  time  T^  in  the  absence  of 

a  driving  source.  The  relaxation  time  Tj  is  a  measure  of  the 

contact  (or  lack  of  contact)  with  the  thermal  bath. 

The  problem  of  spin  relaxation  has  been  studied  by 
36 

Van  Vleck.  His  results  were  not  completely  satisfactory, 

but  the  approach  is  generally  accepted.  The  principles  of  his 
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theory  are  reviewed  here.  Mattuck  and  Strandberg  have 
restudied  the  problem  and  have  arrived  at  substantially  the 
same  results. 

Quantum  mechanically,  thermal  spin  relaxation  is 
treated  as  the  transition  between  spin- states  under  the  excita¬ 
tion  due  to  the  thermal  bath  (in  this  case,  the  lattice).  The 
transition  rate  is  the  reciprocal  of  Tj. 

Spins  can  interact  only  with  magnetic  fields.  The 
lattice,  on  the  other  hand,  is  non-magnetic.  Therefore,  to 
first  order,  there  is  no  thermal  contact  and  T^  is  infinite. 
However,  a  contact  of  the  second  order  does  exist.  In  this 
process,  the  time-varying  electric  field,  due  to  the  thermal 
vibration  of  the  lattice,  interacts  with  the  electrons  in  the 
paramagnetic  ion.  This  affects  the  orbital  moment  of  the  ion 
which,  in  turn,  interacts  with  the  spin  to  induce  the  required 
transition. 

The  total  Hamiltonian  for  the  spin  and  lattice  in  mutual 
contact  is 

H  =  +  (2.13) 
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Hq  is  the  energy  of  the  orbital  state  in  the  crystalline  field 
without  the  spin-orbit  coupling,  is  the  energy  of  the  lattice 
vibrations,  }{  is  the  Zeeman  energy,  Hen  spin-orbit 

energy,  Hqj^  is  the  orbit-lattice  coupling  energy,  and  Hgg  ip 
the  spin-spin  interaction  which  introduces  line  width.  Hq  + 
is  taken  as  the  unperturbed  Hamiltonian,  and  the  remain¬ 
der  is  the  perturbing  Hamiltonian.  It  acts  on  the  manifold  of 
eigenstates  of  Hq  andH^^- 

The  transition  from  the  spin  state  |  +  >  to  |  -  >  is 
computed  by  using  a  "new"  effective  Hamiltonian. 


-ii-)  QC(ij)M'(j.k)H’(k,i')] 

hv(ij)hv(jk) 


(2.14) 


The  indices,  i,  j,  k  and  i'  designate  compound  states  of 
spin,  orbit  and  lattice,  i  and  i'  designate  the  ground  orbital 
state  with  spin  states  |  ->  and  |  +  >,  respectively. 

The  matrix  element  H'(ij)  is  a  sum 


H'(ij)  =  HgQ(ij)  +  HQL(ij)  +  Hjj{ij)  +  Hgg(ij) 


(2.15) 


When  the  ground  orbital  state  is  a  singlet,  HgQ(ij)  and 
Hoj^(ij)  are  zero  whenever  j  contains  the  orbital  ground  state. 
Therefore,  the  sum  over  j  and  k  are  over  states  containing 
only  excited  states.  When  Eq.  (2.  15)  is  substituted  into  £^q. 
(2.  14)  among  the  products  we  find 


Hso(p  +  n.,q  n.)Hso(q-ni,  r-n, )Ho Jr-n.,  p-n.+l) 
[E(p)-E(q)][E(q)-E(r)] 


(2.16) 


and  all  the  permutations  thereof.  Such  terms  represent  third 
order  transitions  from  the  state  |  +  >  to  |  -  > .  At  the  same 
time,  the  lattice  state  changes  from  m  to  m +1.  The  spin  and 
the  lattice  exchange  energy.  Second  order  terms  cancel  out 
almost  completely. 
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The  transition  probability  is  proportional  to  the  square 
of  seen  from  (2.  16)  that  the  probability  is 

proportional  to  .  X.  is  the  spin  orbit  coupling  and  A  is  the 
energy  separation  between  the  orbital  levels  in  a  crystalline 
field.  These  separations  are  not  equal  by  any  means,  but  for 
a  qualitative  discussion,  they  are  all  equated  to  a  single  A. 

The  factor  ^was  previously  encountered  in  the  deriva¬ 
tion  of  the  spin- Hamiltonian.  It  represents  an  orbital  moment 
induced  by  spin-orbit  coupling  in  an  otherwise  "quenched" 
state.  It  was  responsible  for  zero-field  splitting  and  the 
deviation  of  g  from  2.  In  fact,  the  deviation  of  g  from  2  was 
given  by-^.  Therefore,  the  g-factor  is  a  convenient  index 
by  which  relative  relaxation  times  can  be  estimated.  A  long 
relaxation  time  is  incompatible  with  a  large  deviation  of  g 
from  2. 


The  strength  of  the  crystalline  field  enters  into  T 
in  two  ways.  A  is  proportional  to  the  strength  and  so  is 


1 


Hql-  Since 
to 


1 


or  the  transition  probability  is  proportional 

"OL 
■?"  • 


it  appears  that  a  strong  field  results  in  a  longer  relaxation 
time.  The  quartic  dependence  of  the  transition  probability 

on  \  indicates  that  a  very  large  variation  of  T,  is  expected  in 

2+  -1  ^  2+ 
the  iron  group.  V  has  a  \  of  55  cm^  while  Cu  has 

-852cm'f  This  would  indicate  a  ratio  of  relaxation  times 

of  the  order  of  (15)^  50,  000. 
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CHAPTER  III 


THE  RELATIONSHIP  BETWEEN  THE  MATERIAL 
PARAMETERS  AND  MASER  PERFORMANCE 

III.  I  B1  och '  8  Equations 

A  brief  review  is  made  of  the  two-level  maser  to  show 
the  relation  between  parameters  of  the  maser  material  and  the 
maser  performance.  This  relation  is  valid  also  for  three-level 
masers  with  only  a  simple  correction. 

We  have  a  paramagnetic  solid  containing  N  particles  per 
unit  volume,  each  with  a  spin  of  1/2.  In  the  presence  of  a  mag¬ 
netic  field,  the  spin  state  m^  =  1/2  will  have  an  energy  2pH^ 
higher  than  the  state  m^  =  -1/2.  Here  ^  is  the  Bohr  magneton, 
and  H^  is  the  dc  magnetic  field.  Boltzmann's  distribution 
gives  a  population  of  the  lower  state  which  is  larger  than  the 
population  of  the  higher  state  by 


_  A  IT  IKT 

=  ^o  •~-2ffirgT 

where  K  is  Boltzmann's  constant,  and  T  is  the  temperature  in 
Kelvin.  If  AE  <<  KT,  AN?  ^o  the  result 

a  net  induced  magnetic  moment  which  we  shall  call 


Air 

P^  =  PAN=PN^^ 


The  net  angular  momentum  in  units  of'fi  is  =  M^.  The 
analysis  of  a  maser  can  be  made  while  restricting  our¬ 
selves  to  Mq  because  the  remainder  of  the  spins  will  always 
cancel  itself  out. 

The  behavior  of  the  net  moment  M  under  a  time- 
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varying  magnetic  field  excitation  was  analyzed  by  Bloch 
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in  connection  with  nuclear  magnetic  resonance.  The  treatment 
was  made  by  the  now  famous  Bloch  equations.  These  are  classi¬ 
cal  equations.  However,  they  give  results  which  are  in  perfect 
agreement  with  experiments  and  with  quantum  mechanical 
analysis  for  the  simple  cases  where  quantum  mechanical  analysis 
is  possible.  In  addition,  relaxation  effects  can  be  included  in  a 
simple  way;  whereas,  in  a  quantum  mechanical  treatment  it  is 
extremely  complex. 

The  equations  are  basically  the  equations  of  the  spinning 

top. 

=  T  (torque)  =  yTKxE  (3.  2) 

Y  is  the  ratio  of  magnetic  moment  to  the  mechanical  angular 
momentum.  (For  an  electron  this  is  =  .  927x10^.  )  Here,  H 
is  the  total  magnetic  field.  Equation  (3.  2)  can  be  written  for  a 
rectangular  coordinate  system  rotating  about  the  z-axis  with 
an  angular  velocity  ui. 


^  =  -wXS4  +  yKixH  =  i;4X(YR  +  w) 
where  H  is  now 


(3.  3) 


__2  -  ”  f 

a^[H^  +  H^  ^  cos«t]+a^[ — j^cos2wt+ — +a  [ — |^‘sin2ut] 


The  three  components  of  this  equation  are 


dM 


/IT*- 


jr-  —  =  M  [y(H  +H*'  f  cos  wt)  +  w]+M - sin2ut 

dt  y*’  '  o  r.  f.  '  z  2 


dM  H  ,  hi;  2  2 

Y( - ■*’~2 -  ^COSut)+w] 


dM 

W 


yH 


f. 


sin2wt  -  M 


+  H  Coe  2ut)] 
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As  it  stands,  these  equations  are  very  difficult  to  solve. 
A  very  large  simplification  is  made  by  ignoring  the  terms  which 
are  time  varying.  This  would  be  justified  if  the  resulting 
solution  turns  out  to  be  a  slowly  varying  function  compared  to 
the  terms  which  were  ignored. 

The  simplified  equations  are 

dM 


dM 


.1 


H 


r.f. 


dM 


H 


r.  f.. 


The  solutions  are 


(3.4) 


^O^^r.  f. 

M  =  - - - sin  Q  t 

^  2[ar^^ 


M  yH  ,  yH„  +  «  1/0 

o '  r.  f.  '  *  o  '  „ _ ^1/2  ,.v 

M  s - - (1-cos  Q  t) 


(3.5) 


M  1/7  (\H  +w)‘ 

“z'-g-  3- - ">»a  tt— g — 


where  Q  = 


+  (YHp+  w) 


In  the  special  case  of  yH^  +  w  s  0, 


=  0  and 


s  M  cos 
o 


(3.  6) 


This  shows  that  M^  rotates  about  the  x-axis  with  an  angular 
velocity 
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(3.  7) 


w'  = 


vHr. 


f. 


7 


Noting  that  only  one  circularly  polarized  component  of  re¬ 
mains  in  the  solution,  we  can  call  it  H'  to  simplify  the  expressions. 


u'  =  yH'  (3.  8) 

This  is  just  a  free  precession  of  about  the  rotating  field 
H'.  This  condition  is  referred  to  as  magnetic  resonance. 

The  equation  and  solution  just  given  do  not  include 
relaxation  effects.  Bloch  took  relaxation  into  account  by 
inserting  two  damping  times  T^  and  in  the  following 
manner. 

dM  M 

7r^=My(vHo  +  “)-Tf 

dM  M 

gpX  =  M^- yH'-M^(yH^  +  «)-^  (3.9) 

2 


(M  -M  ) 
=  -M  vH'  -  -t'  ° 

y  ^1 


H'  is  the  circularly  polarized  component  of  ^rotating  at  an 
angular  frequency  +<•>  about  the  z-axis. 

These  relaxation  times  are  necessary  for  the  following 
reasons.  First,  there  is  a  continuous  drainage  of  energy  from 
the  spin  system  by  whatever  mechanisms  which  exist  that  are 
trying  to  bring  the  spins  into  thermal  equilibrium.  The  energy 
of  the  spin  system  is  proportional  to  -M  when  the  d.  c.  mag- 
netic  field  is  in  the  z-direction.  Therefore,  a  time  constant 
Tj^  is  inserted  to  point  towards  its  thermal  equilibrium 
value  M^.  Secondly,  the  net  spin  moment  M^  is  not  a  single 
moment  but  is  composed  of  many  spins  oriented  in  the  same 
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direction.  This  orientation  is  imposed  by  the  Boltzmann's 

distribution  when  the  M  is  in  the  z -direction.  But  when  M 

o  o 

is  in  the  x-y  plane,  it  should  relax  towards  the  thermal  equi¬ 
librium  value  for  this  direction  which  is  zero.  In  addition, 
the  dipole-dipole  interaction  between  the  component  spins  try 
to  destroy  M^.  Thus,  for  transverse  components  of  M, 

Bloch  added  a  damping  term 

-  where  <  T^. 


The  solution  to  Bloch's  equation  with  relaxation  is  very 
involved  for  transient  cases.  The  steady  state  solution 
is  very  easy  to  obtain.  It  is 
[(yH  + 

^2  1 

T 


M 


1 


- T ^  T —  - - 7 

^1  ^2  T^  tI 


(3.10) 


M  = 

y 


yH'M^(yH^  +  «) 


TJyH  +w)^+-i2-  +  y^H'^  ] 


yH'M  (yH  +w) 

[(yH  +«)^  +  y‘^H'^  +  -^] 

^2 

at  resonance  yH  +«  =  0.  Then 
*  o 

M  =  0 

X 
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M  = 

y 


YH'M^T2 

- 2 — 2 - 

1  +  yV‘‘TjT2 


(3.  11) 


M  = 
z 


1  +  Y  H'^TjT2 


If  T,  were  infinite,  M  =  M  =  M  =0.  When  T,  is  finite, 

as  it  always  is,  and  become  non-zero.  Physically,  it  is 

easy  to  see  why  must  not  be  zero.  A  finite  means  that 

losses  are  present.  It  indicates  the  presence  of  a  damping  force 

on  SJ.  In  a  steady- state  condition,  an  equal  and  opposite  force 

must  exist  which  balances  out  the  damping  force.  This  is 

provided  by  nxff  =  H  ‘  M  .  For  any  transfer  of  energy  to  take 

X  y 

place  in  the  total  spin  system,  some  transverse  M  must  be 
preserved  even  when  T,  is  present.  The  non-zero  value  of  M 
represents  the  extent  to  which  the  damping  forces  towards 
its  thermal  equilibrium  value. 

The  amount  of  energy  absorbed  by  the  spin  system  from 
the  resonant  magnetic  field  can  be  computed  by  simply  equating 
it  to  the  loss  in  the  system. 

M  -M 
z  o 

P  =  -  ( - iji - )2pH^  =  rate  of  energy  loss  due  to  Tj. 


[  (yH„  +  u))^  +-L-  +  ] 


(3.12) 


The  results  of  Bloch's  equations  can  also  be  expressed 
in  terms  of  susceptibilities. 
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X 


-  JL  - 

^  H  ""  H 

M 


X  = 

DC  ^ 


H 
M  ‘ 


r.f. 


r.f. 


=  2p 


H 


=  0 


=  2p 


M 

3 

H 


Y  *  2P 


(3.13) 


2  1  2  2  ^1 
(yH^+c)  +-  +Y  H-  — 

T2  2 


When  the  standard  phaser  notation  is  used  for  the  time- 

varying  quantities,  the  susceptibility  can  be  represented  as 

a  complex  quantity.  x  =  (x'+jx"),  where  x"  is  a  negative  be- 

1  2 

cause  Y  is  negative.  In  terms  of  the  susceptibility,  P=-jwx"H'  . 

The  solutions  to  Bloch's  equation  have  certain  weaknesses. 

One  is  that  the  equations  do  not  allow  for  the  maintenance  of  a 

38 

transverse  moment.  This  was  pointed  out  by  Redfield.  The 

preservation  of  transverse  moments  for  times  up  to  T.  has 
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been  demonstrated  experimentally.  This  is  of  some 
consequence  in  transient  behavior  when  large  differences 
exist  between  T^  and  but  cause  no  difficulty  in  steady- state 
electron  resonance.  Another  difficulty  is  that  the  power  ab¬ 
sorption  pattern,  as  a  function  of  frequency  which  is  referred 
to  as  "line  shape,  "  is  always  Lorentzian  with  the  Bloch 
equations.  Experimentally,  it  is  not  always  Lorentzian. 

The  value  of  x"  obtained  by  Bloch's  equation  is  inde¬ 
pendent  of  the  magnitude  of  the  individual  moments  which 
make  up  the  total  moment.  This  is  an  e3cpected  weakness  of 
a  classical  treatment.  The  proper  dependence  can  be  found 
by  a  simple  quantum -mechanical  treatment  omitting  explicit 
relaxation  effects. 

From  the  Bloch's  equation  we  had 

2  1  2 
P  =  -2w"H'  =  -  Ti-wx"H 

2  X 

where  H'  was  one  circularly  polarized  component  of  H^.  Therefore, 
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=  2H'.  We  can  consider  this  to  be  the  definition  x".  The  same 
power  P  can  be  written 


P  =  (number  of  transitions)  (energy  of  each  transition) 
x(rate  of  each  transition)  =  AN-  ’  A 


where  A  is  the  transition  probability  of  a  single  ion.  The  quan¬ 
tum  mechanical  transition  probability  is  given  by 


A  =  -^p(k)|  Hj^l  ^ 

where  is  the  matrix  element  of  the  perturbing  Hamiltonian 
which  induces  the  transition  between  states  i|i  and  ,  and 
p(k)  is  the  denstiy  of  states  3=-  .  If  the  density  is  given  as 
^  .  then  • 


A  =  — ^  p(u)  I  H 


kn  I 


For  a  perturbing  Hamiltonian  g^H  *  S 


If  the  spin  states  are  pure,  only  one  component  of  H^,  namely 
H'  contributes  to  the  transition,  giving 


I  ”knl  ^  <k|  S+  I  n>  I  ^ 

p((d)  can  be  interpreted  as  the  line  broadening.  We  can  identi¬ 
fy  this  with 


_ 1 

wTjl 

"^2 


y"h'" 
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which  we  found  in  the  solution  to  Bloch's  equations.  If  A 
represents  the  transition  probability  of  one  ion,  then  p(w) 
must  be  normalized  so  that 

CO 

^  p(u)du  =  1. 

-00 

Integration  of  the  line-shape  factor 


1 

- iji - 

'^2  ^ 
shows  it  to  be  ir.  Therefore, 


p(w)  = 


1 


Substitution  into  P  gives 
Zm  1 


P  =  AN-  AE 


g^P^H'^|<Jc  I  S+  I  n>  I  ^ 


^  *  T2[{YH„-tu>)^+-i,-tYV^] 


2  2 

x"  =  ^N£JL- 


<  k  S+  n  > 


^  T  [(yH  +«)^+-^+yVy^] 

t‘  ^2 


(3.  14) 


The  line  width  of  the  transition  is  given  by  the  factor 
p(u>).  If  the  applied  r.  f.  magnetic  field  is  sufficiently  weak  so 
that  the  last  term  in  the  denominator  can  be  neglected,  the 
half- power  point  in  the  paramagnetic  absorption  is  given  by 
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(yHo  +  “)  =  4^ 

When  the  applied  r.  f.  magnetic  field  is  large  enough  so  the 
last  term  is  not  negligible,  the  half- power  point  is  increased 
to 

(yH^+w) 

The  maximum  absorption  which  occurs  at  resonance  is  also 
reduced  because  of  the  increase  in  the  denominator.  This 
condition  is  called  saturation. 


lU.  2  The  Reflection  Cavity  Maser 


A  maser  can  be  constructed  by  placing  a  suitable 
paramagnetic  material  within  a  reflection  cavity.  The  equi¬ 
valent  circuit  of  such  a  device  is  shown  in  Fig.  3. 


Figure  3 


In  this  circuit,  L  and  C  represent  the  cavity.  R^  represents 
the  resistive  and  dielectric  losses  in  the  cavity  walls  and  in 
the  maser  material.  R^  is  the  absorption  due  to  the  imaginary 
part  of  the  susceptibility,  x"  .  The  relation  between  x"  and  R^ 
can  be  shown  as  follows.  The  power  absorbed  due  to  x"  is 
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The  Q  of  a  cavity  due  to  losses  arising  from  x"  alone  is 


Q  = 


cjX  (energy  stored) 
i’ower  dissipated  by  x" 


-  /  H  dv 

-«x"/H^  dv 


x"  n  (3,  15) 

where  r|  = 

fH^dv 

We  recall  that  H'  is  one  circularly  polarized  component  of 
;  whereas,  H  is  the  total  magnetic  field  intensity.  Thus, 
is  always  less  than  unity.  If  is  the  Q  due  to  alone, 
the  effective  resistance  representing  x"  is 


= 

m 


R,Q 


R^.  2 


(3.16) 


The  total  conductance  of  the  cavity  at  resonance  is 


The  cavity  is  connected  to  a  source  and  load  through 
a  waveguide.  This  coupling  can  be  represented  as  a  trans¬ 
former  as  in  Fig.  4.  The  source  (load)  can  be  transformed 
to  the  secondary  side  to  give  a  new  equivalent  circuit  as  in 
Fig.  5. 

In  this  maser  the  source  and  the  load  are  the  same. 

It  can  be  shown  that  such  a  maser  suffers  from  the  fault  that 
a  large  portion  of  the  source  power  is  dissipated  in  the  source 
(load)  impedance  without  being  amplified  by  the  maser. 
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This  weakness  is  avoided  by  an  insertion  of  a  "circulator" 
between  the  maser  cavity  and  the  load  and  source.  The  circulator  is 
a  three -terminal -pair  device  in  which  energy  flows  only  in  one 
sequence:  a-»  b — ►c.  A  maser  using  the  circulator  is  shown  in  Fig.  6. 

Now  all  the  source  power  is  incident  on  the  maser.  All  the 
output  (reflection)  from  the  maser  is  delivered  to  the  load. 

The  power  delivered  to  the  load  is  obtained  by  simply 
taking  the  reflected  power  from  the  cavity  as  computed  from  trans¬ 
mission  line  equations. 


Figure  4 


Figure  5 
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(3.17) 


z;  -z 

L  o 

Y  -y; 
o  L 

z,  +z 

Y  +y; 

L  o 

o  L 

p.  • 

(3.18) 


For  our  maser  we  have 
1 


Y '  = 
L 


1/R'  -1/R^  -1/R 
o _ m  c 

1/R'  +1/R_  +  1/  R^ 
o  m  c 


(3.19) 


The  condition  for  amplification  is  that 


oo  >  r  >  1 

This  is  achieved  when 


becomes  negative  so  that  the 


numerator  of  Eq.  (3.  18)  is  larger  than  the  denominator.  When 
—  +  j}  is  sufficiently  negative  to  make  the  denominator  vanish, 
the  maser  becomes  an  oscillator. 

The  gain  of  a  maser  was  given  as 


R'  R  R 


R'  is  the  parallel  resistance  of  R  and  R  .  To  obtain  gain,  the 
c  cm  o  » 

following  conditions  are  necessary  or  desirable; 


(1)  <  0 

m 

(2)  I  R  I  <  R 
'  '  I  m '  c 

1  1 _ 1 

^  ^ '  Tr^  r. 


be  as  small  a  positive  quantity 
as  possible 
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It  is  clear  that  to  obtain  amplification  with  a  given  and  R^,  we 
need  a  negative  R^  with  1  sufficiently  small.  This,  in  turn, 

demands  that  x"  be  sufficiently  large,  x",  in  turn,  was  given  by 

_2q  2  2 

(3.  14)  at  resonance:  —  ANp.^2T2'  ^  population  difference 

between  states  1  and  2,  and  matrix  element  of  the  spin 

moment  between  states  1  and  2.  This  shows  that  gain  is  easier  to 
obtain  when 


(1)  The  quantity  AN  is  large.  In  other  words,  we  have 
a  large  inversion. 

(2)  large. 

(3)  T2  is  long. 


Going  back  to  the  expression  for  gain,  we  see  that  another 
way  to  obtain  higher  gain  when  the  conditions  for  amplification  are 
already  satisfied,  is  to  increase  R^.  For  a  given  source  and  load 
impedance,  this  indicates -that  the  coupling  between  the  input-output 
waveguide  and  the  maser  cavity  should  be  made  small.  While  this 
is  certainly  true,  it  restricts  the  bandwidth. 

The  bandwidth  of  the  maser  can  be  limited  by  either  the 
cavity  bandwidth  or  the  paramagnetic  resonance  line  width  of  the 
material.  The  addition  of  the  maser  material  with  negative 
resistance  serves  to  raise  the  Q.  The  expression  for  cavity 
bandwidth  is 


(3.20) 


For  negative  R^ 

According  to  this  expression,  the  cavity  bandwidth  of  the  maser 

cavity  is  a  function  of  both  R'  and  R 
^  o  m 
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The  bandwidth  of  the  empty  cavity,  when  it  is  decoupled 

from  the  load,  is  determined  by  Q^.  This  is  about  5,  000-10,  000. 

The  corresponding  bandwidth  at  10  gc  is  about  2  me.  The  actual 

cavity  bandwidth  for  a  maser  can  be  considerably  narrower  or 

broader.  The  presence  of  broadens  the  bandwidth  while  the 

presence  of  a  negative  reduces  the  bandwidth. 

When  the  bandwidth  of  the  maser  is  limited  by  the  cavity, 
40 

Terhune  shows  that,  although  the  gain- bandwidth  product  is  not 

1/2 

constant,  the  product  (G  -1).  B  is  constant  for  a  given  R  and  R  . 


-  1 


+  ^)=.2f(  > 

m  ^m  c 


[Li/Ql-i/3^+i/3cJ  J 


For  a  negative  R  ,  it  is  2f 
m 


If  the  gain  is  increased  by  increasing  R^,  the  bandwidth  is 
correspondingly  decreased.  The  gain- bandwidth  product  is  increased 
by  decreasing  |  R^  |  . 

If  we  attempt  to  obtain  a  very  broad  bandwidth  for  a  given 

R  ,  R'  must  be  made  low.  However,  a  limit  to  this  approach  is 
m  o 

imposed  by  the  bandwidth  of  the  magnetic  resonance.  The  solution 
for  x"  was  Eq.  (3. 18). 

_ 1 _ 

T2[  ( +  -^  +  ^7 

■^2  ^ 


In  this  expression  resonance  is  given  when  is  set  to  0.  That  is. 


At  a  frequency  w'  i  w. 
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=  Aw,  where  Aw=w'-w 


1 


Tjl  Aw^  +  -ij-  +  ] 

tI  ^ 


2  2  1 

<y  H'  -s —  is  very  small  at  power  levels  applicable  to  normal  maeer 
operation. 

The  magnitude  of  x"  is  -i  of  the  maximum  value  when  A= 


Although  the  line  width  of  the  resonance  and  the  bandwidth  of  the 


reasonable  measure  of  the  bandwidth  of  the  amplifier  when  it  is 


limited  by  the  line  width  of  the  material.  This  is  certainly  the  case 
for  non-resonant  masers  such  as  traveling-wave  masers  which 
we  have  not  discussed.  For  very  wide-band  amplifiers,  T2  must 
be  short. 


Up  to  this  point,  the  relaxation  time  T^  has  not  been  mentioned 
explicitly  in  the  maser  properties.  Yet,  it  is  probably  the  most 
important  parameter  for  maser  materials.  T^  enters  the  problem 
through  the  factor  (N2-N^)  in  the  expression  for  x",  which  was 

already  seen  to  be  the  determining  factor  in  the  gain  of  a  maser. 

4 

Bloembergen's  three-level  scheme,  given  in  the  intro¬ 
duction,  is  to  apply  a  sufficiently  strong  resonant  excitation  be¬ 
tween  the  lowest  and  highest  levels  of  three  energy  levels  so  that 
their  populations,  and  N^,  are  almost  equal.  Thi$  excitation 
is  in  competition  with  T^,  which  represents  a  process  that  is 
trying  to  maintain  thermal  equilibrium. 

The  excitation  power  required  to  make  N2=N^  can  be 
computed.  At  thermal  equilibrium. 


N  AE 

N3  ■  ^1  =  • 

3|C 

When  this  term  is  too  large  to  be  ignored,  it  will  reduce  x", 
and  thereby,  the  gain.  The  maser  is  then  in  nonlinear  operation. 
It  is  then  said  to  be  saturated. 
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When  NjsNj,  the  deviation  of  AN  from  thermal  equilibrium  is 
just  AN.  The  rate  of  energy  loss  through  Tj  is  then 

P=  ae 

2(KT)T, 

A  20 

It  is  instructive  to  estimate  this  quantity.  Suppose  N^  =  10 

AE  =  hx3xl0^°  =  1.  98x10"^^  ergs 

KT  =  5.  35x10“^^  ergs  at  4.  2  °K 
T,=  lO”^ 

Then, 

P  =  335  watts. 

This  is  the  power  which  must  be  made  up  by  the  exciting  field. 

It  must  be  applied  at  a  very  low  impedance  (very  strong  mag¬ 
netic  field)  for  saturation.  This  can  be  seen  from  the  saturation 
condition  in  the  solution  to  Bloch's  equation.  In  the  actual  maser, 
the  cavity  must  be  optimized  for  the  signal  frequency.  We  are  not 
free  to  optimize  the  cavity  for  pumping.  This  factor  and  the 
additional  losses  in  the  cavity  require  that  the  pump  provide 
several  times  as  much  power  as  that  just  estimated.  It  is 
evident  that  the  requirement  is  prohibitive  at  T^  =  10~^.  The 
effect  of  heat  produced  by  such  high  power  dissipation  must  also 
be  considered. 

To  sum  up,  a  short  T^  demands  that  the  pump  power  be 
prohibitively  high  for  maser  operation.  Such  a  high  pump  power 
raises  the  temperature  of  the  maser  even  in  a  thermal  bath. 
Because  T^  decreases  rapidly  with  termperature,  the  pumping 
becomes  inadequate  and  maser  action  cannot  be  obtained. 

Finally,  we  consider  the  question  of  how  the  three 
required  energy  levels  are  to  be  obtained. 
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Confining  ourselves  to  Zeeman  energy  levels,  we  re¬ 
quire  paramagnetic  ions  with  an  angular  momentum  of  at  least 
1  to  give  the  required  three  states  (3=2J+1).  In  the  discussion  of 
the  relaxation  time  T^,  it  was  shown  that  the  presence  of  an 
orbital  moment  leads  to  a  short  relaxation  time.  Therefore, 
we  require  the  orbital  moment  to  be  as  close  to  zero  as 
possible. 

If  the  orbital  moment  were  strictly  zero,  we  are  left 
with  a  pure  spin  state  with  (2  S+1)  -  fold  degeneracy.  If  S  >1, 
the  application  of  a  magnetic  field  will  yield  the  required 
number  of  levels  as  in  Fig.  7.  Unfortunately,  this  scheme 
will  not  work  because  the  spin  states  are  pure.  The  only 
transitions  possible  through  a  time -varying  magnetic  field 
are  those  which  conn^ect  states  whose  values  of  M  differ  by  1. 
All  of  these  transitions  are  resonant  at  the  same  frequency. 

It  is  not  possible  to  pump  between  two  states  which  are  not 
adjacent  to  each  other. 

The  situation  is  rescued  by  a  small  amount  of  orbi¬ 
tal  moment.  It  was  shew  n  that  a  small  amount  of  orbital 
moment  in  the  presence  of  a  sufficiently  anisotropic  crystal¬ 
line  electrostatic  field  will  partly  remove  the  spin  degeneracy 
without  the  application  of  a  magnetic  field.  This  was  called 
zero-field  splitting.  When  the  magnetic  field  is  added,  the 
Zeeman  energy  is  added  to  this  initial  splitting  to  completely 
remove  the  degeneracy.  The  total  result  is  that 

1)  The  resonant  frequencies  for  different  allowed 
transition  are  now  different. 

2)  The  spin  states  can  become  mixed  so  that  many 
mote  transitions  become  allowed.  Pumping  is 
now  possible  between  a  number  of  levels. 
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Energy  levels  and  allowed  transitions 
pure  spin  state.  Example:  S=  5/2. 


Figure  7 


The  Zeeman  energy  is  proportional  to  the  magnetic 
field.  For  a  pure  spin  state  the  relation  is 

AE  =  2p  H 

for  allowed  transitions.  For  a  signal  frequency  of  10  kmc,  the 
field  required  is  3.  5  kilo-gauss.  While  this  is  not  an  unreasonable 
value,  the  magnetic  field  requirement  becomes  prohibitive  for 
very  much  higher  frequencies. 

When  zero-field  splitting  exists,  the  incremental  ratio 
of  frequency  to  magnetic  field  can  be  improved  in  some  cases. 

But,  the  biggest  advantage  of  the  zero-field  splitting  is  that  it 
simply  substitutes  for  Zeeman  splitting.  If  S  is  the  magnitude 
of  zero-field  splitting  given  in  frequency,  it  reduces  the  re- 

g 

quired  magnetic  field  by  "g  f qS  ^^om  the  value  required  for 

pure  spin  at  the  same  frequency  of  resonance.  In  fact,  if  the 

proper  zero-field  splitting  were  avaliable,  no  magnetic  field 

would  be  required.  Masers  based  on  this  scheme  have  been 

41 

discussed  by  Bogle  and  Symons.  It  has  been  successfully 

42  3+ 

operated  by  Kornienko  using  Fe  impurity  in  Al^O^. 

As  the  Zeeman  energy  becomes  large  compared  to 
the  zero-field  splitting,  the  energy  levels  and  states  approach 
the  condition  for  pure  spin,  and  the  difficulties  associated 
with  pure  spin  are  encountered  once  again.  For  this  reason, 
the  upper  frequency  limit  for  a  maser  material  is  given 
approximately  by  the  zero-field  splitting.  Thus,  in  order 
to  extend  the  range  of  the  maser  to  higher  frequencies, 
materials  with  higher  zero-field  splittings  must  be  developed. 

III.  3  Summary 

The  conclusions  of  the  preceding  sections  are  summed 
here  with  some  additional  remarks. 
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The  qualities  desired  of  a  good  maser  material  are 

(a)  A  long  relaxation  time, 

This  is  necessary  if  inversion  is  to  be  accomplished 
without  excessive  pump  power  and  consequent  heating  of  the 
maser  material. 

(b)  A  large  transition  probability  between  the  inverted  states 

The  susceptibility  is  proportional  to  the  square  of  this 
quantity. 

(c)  A  large  density  of  paramagnetic  ions 

Inversion  is  proportional  to  the  density,  other  things 
being  equal.  However,  density  should  be  kept  low  enough 
to  keep  magnetic  ions  separated.  Excessively  high  density 
also  leads  to  another  type  of  relaxation  which  we  have  not 
mentioned:  cross  relaxation.  This  can  lead  to  complications 
-  -  some  good,  some  badf  ^ 

(d)  A  long  T* 

The  susceptibility  is  proportional  to  T2>  There  is  an 
intrinsic  part  which  is  just  inversely  proportional  to  the 
density  of  paramagnetic  ions.  In  many  ions,  there  is  an 
additional  component  due  to  the  nuclear  moments  within 
the  ion  and  in  its  neighbors.  This  should  be  reduced  to  a 
minimum.  Should  a  large  line  width  be  required  for  a 
wide-band  maser,  the  line  can  be  broadened  to  the  desired 
value  by  increasing  the  density  of  paramagnetic  ions.  This 
will  also  help  (c). 

(e)  Zero-field  splitting  equal  to  or  larger  than  the  intended 
maser  frequency 

This  is  necessary  to  make  pumping  possible  and  to 

_ reduce  the  amount  of  magnetic  field  required.  Because  the 

See  the  following  page. 
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T^,  which  is  the  characteristic  time  for  destruction  of  the 
magnetic  moment  normal  to  the  dc  magnetic  field,  is  often 
discussed  in  terms  of  two  separate  effects.  One  is  a  destruc¬ 
tive  interference  of  individual  moments  due  to  their  having  un¬ 
equal  precession  frequencies.  This  can  arise  from  slightly 
different  environments  (position  in  a  magnet  gap,  local  distor¬ 
tions  in  a  crystal,  etc.  )  for  different  spins.  A  slight  varia¬ 
tion  in  precession  frequency  among  the  spins  causes  spins, 
which  were  initially  all  aligned,  to  lose  their  alignment  in 
a  short  time.  This  is  called  inhomogeneous  T2.  The  other 
effect  is  the  resonant  interaction  between  spins  of  the  same 
precession  frequency.  The  interaction  removes  the  degeneracy 
and  thus  broadens  the  line  width.  This  type  of  line  width  is 
described  as  "homogeneous  " 

The  transient  behavior  of  a  spin  system  differs  consider¬ 
ably,  depending  on  whether  the  line  width  is  due  to  homo¬ 
geneous  or  inhomogeneous  However,  in  steady  state, 
as  in  a  three-level  maser,  the  distinction  need  not  be  made. 
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need  for  very  high  frequency  masers  is  most  urgent,  high  zero- 
field  splittings  are  desired. 

CHAPTER  ly 

A  SURVEY  OF  POTENTIAL  CRYSTALS  FOR  THREE-LEVEL 

MASERS 

The  qualitites  desired  in  a  crystal  for  three-level  masers 
were  discussed  in  a  previous  chapter.  The  problem  of  finding 
or  making  crystals  with  those  desired  qualities  is  now  consi¬ 
dered.  This  problem  has  been  discussed  previously  by 
44 

Singer.  There,  only  crystals  which  were  available  syntheti¬ 
cally  were  considered.  Here,  consideration  is  extended  to  all 
crystals.  The  problem  of  finding  or  making  crystals  with 
those  desired  qualities  can  be  divided  into  two  parts.  The  first 
part  is  to  choose  the  paramagnetic  ion  suited  to  the  purpose. 

The  second  is  to  choose  a  suitable  diamagnetic,  dielectric 
solid  to  serve  as  a  host  for  the  paramagnetic  ion.  Of  course, 
the  two  parts  are  not  completely  independent.  The  choice  of 
one  is  bound  to  influence  the  other. 

IV .  1  The  Paramagnetic  Ion 

The  conditions  which  govern  the  choice  are 

1)  Long  Tj- 

2)  Long  T2- 

3)  At  least  three  spin  levels. 

Van  Vleck's  theory,  showed  that  spin-orbit  coupling  was 
the  mechanism  which  was  responsible  for  T^.  A  long  Tj^ 
requires  a  small  spin-orbit  coupling  constant  \  and  a  small 
orbital  moment.  States  which  belong  to  terms  of  L=0  satisfy 
this  requirement.  In  addition,  even  where  L  ^  0,  crystalline 
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fields  of  the  proper  symmetry  were  found  to  split  the  states  in 
the  manifold  of  L  in  such  a  way  that  left  a  singlet  state  with  the 
lowest  energy,  fn  such  a  singlet  state,  the  only  orbital  moments 
present  were  those  re-induced  by  spin-orbit  coupling.  This  was 
a  second-order  process  proportional  to  —  ,  where  Aj  is  the 
energy  from  the  ground  state  to  a  higher  orbital  state  j.  There¬ 
fore,  in  such  cases,  in  addition  to  a  small  X.,  large  A.'s  are 
required  for  long  T^.  This  second  condition  dictates  that  the 
ground  orbital  state  be  a  singlet  in  an  octahedral  or  cubic 
syrhmetry.  This  is  because  the  common  environments  of  the  iron 
group  are  usually  very  close  to  one  of  these  two  symmetries. 
The  exact,  actual  symmetry  may  be  lower  because  of  small 
deviations,  but  the  removal  of  degeneracy  by  such  lower  sym¬ 
metries  will  not  result  in  A's  large  engouh  to  produce  good 
quenching. 

An  exception  is  made  for  "non -magnetic"  doublets  dis¬ 
cussed  in  Section  II.  9.  The  component  states  in  such  a  doublet 
have  quenched  orbital  moments.  In  addition,  there  is  no  first 
order  mixing  between  the  components  under  spin-orbit  coupling 
so  that  the  quenching  is  preserved. 

The  strong  dependance  of  T^  on  X  must  be  noted  in  selecting 
ions  for  long  T^. 

With  respect  to  T^,  the  following  combinations  of  ion  and 
symmetry  are  suitable. 


Ion 

Ti3+ 


Cr 

Cr 


3+ 


2+ 


Mn 


2+ 


Fe 

Fe 

Co 


3+ 

2+ 

2+ 


M-2+ 

Ni 


Symmetry 

cubic 

cubic 

octahedral 
octahedral 
octahedral 
any  symmetry 
any  symmetry 
cubic 
cubic 

octahedral 

octahedral 


Orbital  Ground  State 

non- magnetic  doublet 
singlet 
singlet 
singlet 

non-magnetic  doublet 
L  =  0 
L  =  0 

non-magnetic  doublet 
singlet 
singlet 

non-magnetic  doublet 


See  the  following  page. 
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3+  2+ 

Among  the  combinations  listed,  Fe  and  Mn  are 
superior  to  the  others.  Cr^^  in  an  octahedral  field  is  also 
better  than  the  others  because  of  the  large  energy  gap  between 
the  singlet  and  the  next  lowest  orbital  and  a  small  spin- orbit 
coupling.  Trivalent  ions  have  larger  crystal  field  splitting  and 
are,  therefore,  to  be  preferred.  The  non-magnetic  doublets 
are  generally  inferior  to  the  singlets.  Cu^'*’  and  Ni^^  have 
rather  large  values  of  \.  Therefore,  they  are  less  desirable. 

The  contribution  by  the  ion  to  T2  is  the  inhomo¬ 
geneous  broadening  caused  by  the  nuclear  spin.  When  the 
hyperfine  splitting  is  large,  the  lines  outside  the  intended 
bandwidth  of  the  maser  do  not  participate  in  the  maser  action 
and  is  wasted.  The  effective  density  of  the  spins  is  reduced. 

The  choice  of  ions  with  respect  to  T,  falls  on  those  ions  with 

^  2+  3+ 

a  minimum  of  nuclear  spin.  The  desirable  ions  are  Fe  ' 

Ni^”*"  and  Cr^^’  in  the  order  given. 

The  condition  S  >  1  is  also  satisfied  by  Fe' 

Ni^'*‘  and  Cr^'*'' 

Ions  which  satisfy  all  the  conditions  are  Fe^**^  in 
2+  3+  2+ 

any  symmetry,  Cr  ,  Cr  and  Ni  in  octahedral  symmetry 

2+ 

and  Fe  in  cubic  symmetry. 

In  addition  to  the  iron  group  ions,  the  two  S- state 

3+  2+ 

rare  earths,  Gd  and  Eu  ,  display  long  T^*  The  other  rare 

earths  have  short  relaxation  times  due  to  excessive  orbital 
2+ 

moment,  Eu  is  not  suited  for  masers  because  of  its  nuclear 
moment. 

The  useful  ions  listed  in  the  probable  order  of 
Cr^+,  Gd^+,  Ni^+,  Cr^+  and  Fe^^. 

^By  cubic  symmetry,  8  negatively  charged  neighbors  in  a 
cube  Are  designated.By  octahedral  symmetry,  6  nej;atively 
charged  neighbors  in  an  octahedron  are  designated. 


2+,  3+ 
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IV.  2  The  Host  Crystal 

The  desired  qualities  in  a  host  crystal  are 

a)  Absence  of  paramagnetism  for  ferromagnetism  on  the  part 
of  the  pure  host  crystal  itself. 

b)  A  minimum  abundance  and  strength  of  nuclear  moment,  par¬ 
ticularly  in  the  ions  immediately  adjacent  to  the  paramagnetic 
impurity. 

c)  Sufficient  asymmetry  to  produce  a  large  zero-field  splitting. 

d)  Good  stability.  Particularly,  resistance  to  thermal  stress 
and  moisture. 

e)  A  small  unit  cell,  so  that  all  the  paramagnetic  ions  are 
magnetically  equivalent. 

f)  If  possible,  cations  of  the  same  valence  as  the  intended  para¬ 
magnetic  substitute. 

g)  Cations  whose  ionic  radii  are  similar  to  the  radii  of  the 
intended  paramagnetic  substitute. 

h)  If  possible,  sufficiently  low  symmetry  to  remove  orbital 
degeneracies  in  the  ground  state. 

i)  Very  high  resistivity. 

The  condition  on  nuclear  moment  rules  out  a  very  large 
number  of  the  most  commonly  available  crystals.  All  the 
halides  are  ruled  out  because  of  their  nuclear  spin.  Also  ruled 
out  are  all  the  hydrated  crystals  because  of  the  nuclear  spin 
of  hydrogen.  For  the  same  reason,  most  organic  crystals 
are  not  suitable. 

The  crystals  which  may  be  suited  are  the  oxides,  car¬ 
bonates,  sulphates,  tungstates,  chromates,  silicates,  silicide^, 
sulphides  and  carbides  of  metals  which  are  not  in  themselves 
paramagnetic  and  which  have  a  low  abundance  of  nuclear 
moment. 

The  condition  that  large  zero-field  splitting  should  result 
is  not  as  easy  to  apply.  The  theoretical  relation  between  the 
zero-field  splitting  and  the  strength  and  symmetry  of  the  crys¬ 
talline  field  is  well  known,  but  so  far,  the  numerical  values  of 
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the  parameters  which  enter  into  the  expression  are  not  known 
with  sufficient  accuracy  to  allow  even  approximate  calculation. 

We  do  not  know,  however,  that  a  cubic  symmetry  cannot  cause  a 
zero-field  splitting  in  S=1  or  S=3/2,  Therefore,  in  the  case  of  the 
paramagnetic  ions  Cr^^  or  Ni^^,  a  cubic  symmetry  must  be 
avoided.  A  word  of  caution  is  in  order  here.  The  symmetry, 
which  in  involved,  is  the  point  symmetry  of  the  paramagnetic 
ion  site  and  not  of  the  crystal  as  a  whole.  These  two  symmetries 
usually  differ,  the  point  group  being  lower. 

It  was  noted  previously  that  the  positions  of  the  nega¬ 
tively  charged  neighbors  about  the  metallic  cation  are  usually 
very  close  to  octahedral  when  the  radius  of  the  metallic  ion 
is  small  compared  to  the  radius  of  the  negative  neighbors. 

When  the  negative  neighbor  is  oxygen,  this  condition  is 
usually  satisfied.  We  can  speculate  that  this  relation  is  not 
followed  as  closely  in  the  case  where  the  neighboring  oxygens 
belong  to  tightly-bound  complexes  such  as  SiO^,  SO^,  CO^, 

CrO^  or  WO^.  In  such  cases,  the  oxygens  may  not  be  free 
to  pack  themselves  as  closely  and  regularly  about  the  metallic 

ion.  The  regularity  can  also  break  down  when  the  metallic 
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ion  is  comparable  to  or  larger  than  oxygen  in  radius. 

The  unit  cell  should  be  small  in  order  that  the  para¬ 
magnetic  ions  are  not  placed  in  a  number  of  sites  which  do  not 
give  the  same  spectrum.  Such  a  case  reduces  the  effective 
number  of  ions  participating  in  any  one  resonance  except  for 
special  orientations  of  the  magnetic  field.  By  a  large  unit  cell 
we  mean  that  the  unit  cell  contains  a  large  number  of  molecules. 
Magnetically,  many  of  these  may  be  equivalent.  Two  sites 
related  by  inversion,  for  example,  are  always  equivalent 
magnetically  because  of  the  inversion  symmetry  of  all  para¬ 
magnetic  spectra.  The  number  of  nonequivalent  sites  is  also 
reduced  when  the  spin  value  S  is  small  or  the  symmetry  is 
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very  high,  for  we  have  seen  that,  in  these  cases,  details  of  the 
symmetry  fail  to  be  registered  in  the  magnetic  spectra.  If  the 
sites  in  a  unit  cell  are  nonequivalent  only  because  of  these 
details  which  cannot  be  registered  in  the  paramagnetic  spectrum, 
the  sites  become  equivalent  magnetically.  In  many  cases,  even 
when  the  sites  are  not  completely  equivalent,  there  are  ample 
special  orientations  which  make  them  equivalent  to  make  them 
useful  in  maser  work. 

The  requirement  of  proper  valence  is  not  a  stringent  one. 

It  is  true  that  extensive  replacement  of  the  cation  by  single 

impurity  requires  the  valence  state  of  the  replacement  ion  to 

be  the  same  as  the  original  ion.  For  dilute  substitutions  of  less 

than  .1%,  this  is  not  necessary,  particularly  if  the  crystal  is 

not  annealed.  The  excess  or  deficiency  of  charge  can  be 

neutralized  by  a  compensating  impurity  or  by  an  additional 

anion  in  the  nearest  interstitial.  If  the  interstitial  anion  or 

charge  compensating  second  impurity  is  close  enough  to  the 

paramagnetic  ion,  the  symmetry  about  the  ion  is  altered  suf- 

28 

ficiently  to  influence  the  spin  resonance  spectrum.  Tinkham 
has  observed  such  spectra  for  Cr^"*"  in  ZnF2.  Very  often, 
however,  the  compensation  is  sufficiently  far  away  so  that 
no  detectable  change  occurs  in  the  symmetry  about  the 
paramagnetic  ion. 

Until  recently,  it  was  believed  that  substitutions  could 

not  occur  if  the  ionic  radii  of  the  ions  to  be  exchanged  differed 

by  more  than  15%.  However,  the  research  in  maser  materials 

within  the  last  two  years  has  shown  that  substitutions  up  to  the 
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order  of  .  1%  may  be  possible  in  many  cases. 

The  requirement  that  the  resistivity  be  very  high  is 
obvious.  Even  a  small  resistivity  would  greatly  reduce  the 
Q  of  the  maser  cavity.  Furthermore,  at  microwave  fre¬ 
quencies,  conductivity  leads  to  skin  effect  which  removes  the 
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interior  of  the  crystal  from  participation  in  maser  action.  In 
this  respect,  the  sulphides  must  be  examined  carefully.  Some 
of  the  sulphides  are  fairly  good  conductors  or  semiconductors. 
Galena  (PbS)  is  a  good  example. 

Some  of  the  sulphides  have  been  used  as  hosts  for  para¬ 
magnetic  impurities.  Mn^^  in  ZnS  and  Eu^^  in  SrS  are  notable 
examples.  Sulphides  are  not  always  ionic,  and  many  of  them 
may  prove  unsuitable  for  masers.  The  same  situation  exists 
with  carbides  and  silicides. 

Armed  with  the  considerations  stated  above,  we  can 

examine  the  compilation  of  inorganic  compounds  given  in  the 

47 

Handbook  of  Chemistry  and  Physics,  WycTcoff's  Crystal 
Structures  and  the  International  Tables  for  X-Ray  Crystallo- 
graphy.  Two  facts  stand  out  in  such  an  examination. 

a)  The  lack  of  suitable  crystals  which  give  a  cubic 
(8-fold  coordination)  symmetry.  The  CaF2  (fluorite)  structure 
is  the  only  such  structure.  All  but  one  member  of  this  class 
must  be  rejected  for  one  reason  or  another.  CeO^  is  the  only 
possibility.  Unfortunately,  it  is  accurately  cubic.  In  addition, 
its  quadri-valent  cation  makes  it  less  useful.  Three  other 
structures  have  8-fold  coordination,  but  they  are  very  far  from 
cubic.  They  are  the  ZrSiO^  (zirzon)  structure,  CaWO^ 
(scheelite)  structure  and  Zr02-  Of  these,  ZrSiO^  and  CaWO^ 
are  very  interesting  for  other  reasons.  Zr02  is  difficult  to 

use  because  of  the  quadri-valent  cation  and  its  four  nonequi¬ 
valent  sites. 

b)  The  remarkably  small  number  of  suitable  crystals 
with  trivalent  cations.  L.a202  is  practically  the  only  promising 
one  besides  the  alr^^ady  well-known  Al20^  (corundum). 

A  list  of  some  crystals  which  may  be  suitable  is  given 
below.  This  is  by  no  means  complete.  On  the  other  hand,  many 
of  these  crystals  may  be  unsuitable  because  of  some  factor 
which  escaped  the  author's  consideration. 
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BaCOj!  BaCrO^,  BaO,  BaSO^,  BaWO^,  BaSiO^ 

CdC03.  CdSiO^,  CdSO^,  CdWO^ 

OaOO^y  CaSiO^i  CaCrO^ 

Ce02 

PbC03,  PbCrO^,  PbSO^,  PbO,  PbSO^,  PbSi03 
PbWO^,  MgO,  MgC03,  MgWO^,  MgSi03,  Si02 
SrC03,  SrCrO^,  SrO,  SrSi03,  SrSO^,  SrWO^ 
SnO,  Sn02,  SnSO^ 

Ti02 

ZnCrO^,  ZnC03,  ZnO,  ZnSi03,  ZnSO^ 

Zr02,  ZrSiO^ 


The  exact  structures  of  many  of  the  silicates  are  not  yet 
known.  The  metasilicates  of  known  structure  have  very  large  unit 
cells,  CaSi03,  for  example,  has  12  molecules  in  a  unit  cell. 

Perhaps,  the  only  interesting  crystal  in  the  group  is  zircon,  which 
has  four  molecules  in  a  unit  cell.  Fortunately,  all  the  sites  are 
equivalent.  The  local  symmetry  is  42m,  8  oxygens  form  the 
nearest  neighbors. 

The  sulphates  which  are  listed  are  all  of  the  same 
orthorhombic  structure.  There  are  four  molecules  in  a  unit  cell. 
However,  they  are  made  up  of  pairs,  related  by  inversion  so  that 
only  two  sites  are  nonequivalent  magnetically.  The  radii  of  the 
cation  are  very  large.  This  results  in  a  high  coordination  number (12) 
and  low  point  symmetry  (m).  All  the  cations  are  divalent. 

The  carbonates  are  of  two  distinct  structures.  The  smaller 
cations  result  in  the  trigonal  structure  of  CaC03  (calcite).  The 
point  symmetry  at  the  cation  site  is  also  trigonal  (C3^),  The  two 
sites  in  the  unit  cell  are  not  equivalent.  The  structure  of  the  car¬ 
bonates  with  large  cations  is  similar  to  the  sulphates.  The  point 
symmetry  about  the  cation  is  m.  The  cations  are  surrounded  by 
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six  oxygens  in  a  distorted  octahedron. 

The  tungstates  have  already  been  successfully  exploited 
as  hosts  for  paramagnetic  ions.  They  are  of  two  types.  Those 
with  large  cations  form  tetragonal  crystals.  There  are  four 
molecules  in  a  unit  cell,  but  all  are  equivalent.  This  is  a  parti¬ 
cular  virtue  of  this  structure.  Tungstates  with  small  cations  form 
monoclinic  crystals  which  are  isomorphous  (have  the  same 
structure  and  size;  therefore,  interchangeable)  with  most  of  the 
tungstates  of  the  iron  group.  They  have  two  molecules  in  a 
unit  cell.  The  two  cation  sites  are  equivalent.  Both  types  have 
shown  great  promise  in  maser  applications. 

The  chromates  are  divided  into  three  structures.  Those 

with  big  cations  (CaCrO^)  have  the  same  structure  as  ZrSiO^ 

and  are  very  suitable.  The  four  cation  sites  in  a  unit  cell  are  all 

equivalent.  The  local  symmetry  is  42m.  Another  structure  is  that 

of  ZnCrO.  and  CdCrO..  These  crystals  are  orthorhombic.  The 
’  2  ’ 

local  symmetry  is  —  .  The  four  cations  in  a  unit  cell  are  all 
equivalent.  The  coordination  around  the  cation  is  a  distorted 
octahedron.  PbCrO^  is  monoclinic  and  has  two  nonequivalent 
sites. 

The  oxides  are  also  divided  into  a  small  number  of  groups. 
The  oxides  MgO,  CaO,  BaO,  SrO  and  CdO  are  accurately  cubic 
with  exact  octahedral  coordination  about  the  cation. 

A  very  interesting  crystal  is  La^p..  It  has  only  one  mole- 

3+  ^  ^ 

cule  in  a  unit  cell.  The  two  La  in  the  unit  cell  occupy  sites 
which  are  related  by  inversion,  making  the  sites  equivalent 
magnetically.  The  local  symmetry  is  3m.  Furthermore,  it 
is  an  ideal  host  for  trivalent  ions.  Gd^^  in  La^O^  should  be 
an  ideal  combination. 

Ti02  and  Sn02  constitute  another  class.  The  unit  cell 
contains  two  cation  sites  which  are  not  equivalent.  The  local 
symmetry  is  orthorhombic.  The  apparent  success  of  Ti02  as 
a  maser  material  raises  hopes  that  other  crystals  with 
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quadri-valent  cations  may  be  successful.  CeO^  with  the 
accurately  cubic  fluorite  structure  may  be  such  a  crystkl. 

PbO  and SnO have  atetragonal  crystal  with  2  cations  per  unit 
cell.  The  local  symmetry  is  given  as  4mm.  The  two  sites  are 
equivalent  niagnetically.  This  crystal  is  not  strictly  ionic,  and 
its  suitability  is  a  question  mark. 

CONCLUSION 

The  following  combinations  of  crystal  and  paramagnetic 
impurity  seem  to  be  particularly  promising. 

La-O,  It  should  be  an  ideal  host  for  Gd^^.  It  should  also  be 
^  ^  3+ 

very  good  for  Fe  .  However,  the  large  difference 
in  ionic  radii  causes  some  concern.  The  7-fold  co¬ 
ordination  does  not  permit  us  to  say  what  the  orbital 
ground  state  of  Cr^^,  Ni^^,  Cr^^  and  Fe^^  will  be 
with  any  assurance.  They  may  be  favorable.  Rare 
earths  ar6  difficult  to  separate  from  each  other. 
Growing  pure  La^O^  without  other  rare-earth 
impurities  may  be  difficult.  The  line  broadening 
due  to  the  nuclear  moment  of  La  should  not  be  too 
large  because  La  would  not  be  the  nearest  neighbor 
of  paramagnetic  impurities. 

CaCrO.  Fe^^  and  Gd^^  should  be  good.  Again,  the  large 

’  2+ 

difference  in  ionic  radii  of  Ca  and  the  iron  group 

impurities  causes  some  concern.  The  8-fold  coordi¬ 
nation  is  far  from  cubic  so  that  it  is  difficult  to  say 
what  the  orbital  ground  states  for  Cr^^,  Ni^^,  Fe^^ 
and  Cr^^  will  be.  They  may  be  favorable.  Irregular 
coordination  may  result  in  large  aero-field  splitting. 
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ZnCrO . 
4 

CdCrO. 

4 


ZrSiO^ 

CaWO.'l 

4 

BaWO^ 
SrWO^  . 


MgWO^'i 
ZnWO^  I 

CdWO^  J 


In  spite  of  having  four  molecules  in  a  unit  cell,  all  sites 
are  equivalent.  The  point  of  symmetry  is  only  mm.  The 
4-fold  coordination  about  the  Zn,  Cd,  sites  make  them 
particularly  attractive.  In  some  structures  the  positions 
of  Cr^^  and  the  other  metallic  cation  become  interchanged. 
The  situation  is  identical  with  CaCrO.. 

4  3. 

These  crystals  have  already  been  examined  with  Fe 
3+  13 

and  Gd  The  results  have  been  excellent.  The 

success  in  introducing  Fe^^  gives  hope  that  other 
crystals  with  large  cations  will  likewise  accept  Fe^^. 

The  coordination  is  8-fold  but  far  from  cubic.  The 
absence  of  any  report  on  Cr  may  indicate  that  the 
sign  of  the  crystalline  field  puts  the  orbital  triplet 
lowest. 

These  crystals  have  shown  very  good  characteristics 

^  3+  .  12 

with  Fe  impurity. 


The  following  crystals  will  yield  two  spectra  for  ions 
with  S  >  2,  Nevertheless,  they  show  features  which  warrant  investi 
gatibn. 


CaCOj'i 

MgC03 I 
ZnCO^  J 


BaCO^ ' 

SrCO^  ► 

PbCO^ , 

PbCrO. 

4 


These  crystals  should  be  good  hosts  for  Fe^^,  Gd^'*’, 

2+  2+ 

Ni  and  Cr  .  The  coordination  is  close  to  octahedral 
assuring  a  singlet  orbital  state  for  Ni^^  and  Cr^^. 
Structure  is  fairly  similar  to  AI2O3.  The  nuclear  spin 
content  is  exceptionally  low. 

These  crystals  also  feature  low  nuclear  moment  and  an 
octahedral  coordination. 
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BaSO^  > 
SrSO^  I 

PbSO^J 


Ti02 

Sn02 


MgO 

BaO 

CaO 

Ce02 


The  very  high  (12)  andir  regular  coordination  makes  it 
difficult  to  make  any  predictions.  If  or  Gd^^ 

can  be  introduced,  large  zero-field  splitting  may 
result.  It  is  difficult  to  predict  the  degeneracy  or 


near  degeneracy  of  the  orbital  ground  state  for  Cr 

r'  2  + 

Cr 


'  3+ 


Fe^"^  or  Ni^'*'. 


3+ 


Ti02  crystals  with  Fe  have  already  shown  huge 
zero-field  splitting. 


Other  interesting  combinations  are  as  follows: 


The  pure  octahedral  field  can  give  zero-field  splitting 

in  Fe^^  and  Gd^^.  Nevertheless,  the  high  symmetry 

makes  them  difficult  to  use  in  three-level  masers. 

This  is  the  only  crystal  with  8-fold  coordination 

2+ 

cubic  symmetry  about  the  cation.  Fe  should  be 
very  interesting  if  it  can  be  made  to  substitute  the 
quadri-valent  Ce.  Very  large  zero-field  splitting 
should  result. 


Before  concluding  this  section,  two  more  classes  of 
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crystals  deserve  mention.  They  are  the  borates  and  the  phos¬ 
phates.  These  crystals  were  excluded  from  the  preceding  dis¬ 
cussion  because  of  their  abundant  nuclear  spin.  However,  as  in 
AI2O3,  the  nuclear  spins  reside  not  in  the  immediate  neighbors 
but  in  the  next  nearest  neighbors  of  a  paramagnetic  impurity. 
Therefore,  the  contribution  to  the  line  width  is  small.  The  borates 
are  particularly  attractive  because  they  form  stable  crystals, 
and  they  generally  form  crystals  with  trivalent  cations.  There 
is  a  lack  of  literature  on  the  borates.  However,  it  is  found  that 
LaBO^  has  the  same  structure  as  calcite.  It  should  be  an  ideal 
host  for  Fe^^,  Gd^^  and  Cr^^.  The  approximately  octahedral 
arrangement  of  oxygens  around  the  cation  should  give  a  singlet 
orbital  state  as  the  ground  state  for  Cr^^. 
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IV.  3  Availability  of  Crystals 


Most  of  the  crystals  in  the  preceding  section  have  been 
known  as  natural  crystals;  some  as  gems  and  semiprecious 
stones.  Various  impurities  are  present  in  natural  crystals, 
and  paramagnetic  resonance  due  to  such  impurities  have  been 
studied.  Only  in  a  few  fortunate  instances,  however,  do  the 
crystal  and  impurity  form  the  combination  desired  in  a  three - 
level  maser.  Even  when  the  right  combination  is  present,  the 
concentration  may  not  be  optimum;  an  abundance  of  other  un¬ 
desirable  impurities  may  be  present;  or  single  crystals  may 
be  too  small.  A  crystal  must  be  available  synthetically  before 
it  can  be  a  practical  maser  material.  Only  then  can  all  the 
properties  be  controlled  to  specifications.  The  particular 
prominence  of  ruby  (Al^O^iCr)  as  a  three-level  maser 
material  is  largely  due  to  its  availability  in  synthetic  form. 
Ti02  is  also  investigated  because  it  is  available  synthetically. 

Growing  synthetic  crystals  is  not  necessarily  difficult. 

In  a  great  many  cases,  it  is  simply  the  lack  of  sufficient 
scientific  or  commercial  incentive  which  has  resulted  in 
their  neglect.  Perhaps,  in  a  number  of  crystals,  all  that 
is  required  is  to  try.  However,  the  very  great  difficulty 
experienced  in  the  growth  of  synthetic  diamond  and  quartz 
should  restrain  one  from  too  much  optimism.  Generally 
speaking,  crystals  which  can  be  melted  at  reasonable 
temperature  and  pressure,  or  crystals  which  can  be  dis¬ 
solved  in  some  solvent  are  probably  easiest  to  synthesize. 

For  example,  hydrated  crystals  can  be  grown  readily 
from  solution.  The  very  insoluble  BaSO^  is  difficult  to  grow. 
CaWO^  can  be  grown  from  a  melt.  Crystals  such  as  CaCO^ 
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decompose  upon  heating  so  that  growth  from  a  melt  is  not 
possible. 

Among  the  crystals  which  were  classified  as  promising 
or  interesting  in  the  preceding  section,  the  following  have 
been  grown;  MgO,  CaWO^,  MgWO^,  ZnWO^,  CdWO^,  and 
TiO^. 
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